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1.  Introduction 


To  understand  the  deflagration  of  ammonium  perchlorate-hydroxyl-terminated 
polybutadiene  (AP-HTPB)  composite  propellants  at  a  fundamental  level,  US  Army 
Research  Laboratory  (ARL)  researchers  have  built  models  to  simulate  the  process 
(Chen  and  McQuaid  2014;  McQuaid  and  Chen  2014;  Nusca  2014).  Of  particular 
interest  is  the  mechanism  (or  mechanisms)  underlying  the  abrupt  change  in  the 
slope  of  the  burning  rate  (n)  versus  pressure.  This  phenomenon,  which  typically 
occurs  around  3,000  psia  (Atwood  et  al.  1999)  and  is  referred  to  as  an  “exponent 
break,”  limits  the  performance  potential  of  this  propellant  class  (Atwood  et  al. 

2013) . 

At  the  start  of  our  efforts  (in  2012),  we  found  that  the  best  available  finite-rate 
chemical  kinetics  mechanism  for  modeling  the  gas-phase  reaction  chemistry 
associated  with  AP’s  combustion  had  a  number  of  questionably  parameterized  rate 
expressions  that  sensitivity  analyses  indicated  were  important.  Addressing  this 
issue  through  the  application  of  computational-chemistry-based  methods  and 
thermochemical  kinetic  theories,  we  developed  and  validated  a  gas-phase  finite- 
rate  chemical  kinetics  mechanism  for  modeling  the  combustion  of  pure  AP  (Chen 
and  McQuaid  2014).  That  mechanism  has  since  been  updated  and  integrated  with 
a  previously  built  chemical  kinetics  mechanism  for  the  combustion  of  HTPB 
(R45M  polymer)  in  an  NOx-rich  environment  (Chen  and  McQuaid  2009;  Chen  and 
McQuaid  2010),  producing  a  mechanism  for  modeling  the  combustion  chemistry 
associated  with  AP-HTPB  deflagration.  It  has  (at  present)  2,627  elementary 
reactions  and  involves  637  species.  Some  details  of  the  estimation  of  reaction-rate 
expressions  and  thermochemical  property  estimates  that  were  developed 
specifically  for  this  mechanism  have  been  presented  (Chen  and  McQuaid  2015). 

Because  simulation  run  times  would  be  prohibitively  long  and  memory 
requirements  excessive,  it  is  not  practical  to  employ  the  full  AP-HTPB  mechanism 
as  a  submodel  in  computational  fluid  dynamics  (CFD)  frameworks  that  our  group 
employs  to  simulate  deflagration  phenomena  (McQuaid  and  Chen  2014;  Nusca 

2014) .  Therefore,  the  trial  mechanism  method  (TMM)  invented  by  Kotlar  (2010) 
and  further  developed  by  McQuaid  (2013)  was  employed  to  produce  (reduced) 
skeletal  sets.  With  an  objective  of  producing  a  mechanism  with  100  or  fewer 
reactions  and  involving  100  or  fewer  species  to  study  exponent  break  phenomena, 
screening  criteria  were  established  to  indicate  a  candidate’s  ability  to  mimic  the  full 
mechanism  in  models  representing  the  deflagration  of  AP-HTPB  mixtures  having 
from  80-  to  90-wt%  AP.  The  pressure  range  of  primary  interest  was  1,000-10,000 
psia. 
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As  in  previous  mechanism  reduction  efforts  based  on  the  TMM  (Kotlar  2010; 
McQuaid  2013;  McQuaid  and  Chen  2014),  results  of  homogeneous  reactor  (HR) 
model  simulations  produced  with  trial  mechanisms  were  employed  as  bases  for 
eliminating  reactions.  However,  the  implementation  of  the  method  was  not  as 
straightforward  as  it  was  in  the  past.  Because  AP-HTPB  propellants  have 
composite/nonhomogeneous  structure,  there  are  regions  near  the  condensed-phase- 
gas-phase  interface  that  are  very  rich  in  the  daughters  of  HTPB’s  pyrolysis  and 
others  that  are  rich  in  the  daughters  of  AP’s  decomposition  (Surzhikov  and  Krier 
2003;  Gross  et  al.  2013).  As  such,  there  was  a  desire  to  screen  trial  mechanisms  on 
the  basis  of  simulations  with  starting  chemical  compositions  that  ranged  from  the 
nascent  product  distributions  for  the  condensed-phase  to  gas-phase  conversion  of 
(pure)  AP  to  the  nascent  product  distribution  for  the  condensed-phase  to  gas-phase 
conversion  of  (pure)  HTPB  (R45M  polymer)  (Surzhikov  and  Krier  2003;  Gross 
et  al.  2013).  Vetting  trial  mechanisms  at  a  number  of  different  pressures  over  the 
range  of  interest  was  also  desired.  However,  because  the  run  times  of  relevant 
simulations  were  relatively  long,  the  number  of  simulations  that  could  be  employed 
for  screening  was  extremely  limited.  Since  that  limited  the  parameter  space  within 
which  candidates  produced  by  the  TMM  could  be  expected  to  be  valid,  it  prompted 
us  to  make  some  assumptions  about  the  relative  importance  of  various  subprocesses 
and  to  apply  the  method  in  stages.  Specifics  of  the  implementation,  including  some 
deviations  from  past  practice,  are  discussed. 

Postreduction  analyses  were  conducted  to  identify  candidates  with  the  best 
potential  to  mimic  the  full  mechanism  in  intended  applications.  Of  the  candidates 
with  the  potential  to  be  a  submodel  for  the  ARL  CFD  laminate  framework  (Nusca 
2014),  the  one  considered  to  be  the  best  compromise  between  size  and  ability  to 
mimic  the  full  mechanism  within  the  delineated  parameter  space  had  100  reactions 
and  involved  81  species.  Results  from  simulations  that  were  the  basis  for 
mechanism  reduction  and  candidate  selection  are  presented.  Burning  rate  estimates 
for  homogeneous  AP  and  AP-HTPB  systems  produced  with  a  1 -dimensional  (1-D), 
2-phase,  premixed  laminar  flame  (CYCLOPS)  model  provide  additional  evidence 
for  the  mechanism’s  validity  for  intended  applications. 

2.  Implementation  of  the  Trial  Mechanism  Method 

More  complete  descriptions  of  the  TMM  are  given  elsewhere  (Kotlar  2010; 
McQuaid  2013).  Steps  include  the  following: 

.  Randomly  ordering  the  full  mechanism’s  reactions/rate  expressions. 

.  Sequentially  eliminating  single  reaction  rate  expressions  from  it  on  a  trial 
basis. 
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Rerunning  constant-pressure  (CONP)  HR  simulations  with  the  (trial) 
mechanism  created  by  the  elimination. 


.  Permanently  eliminating  the  reaction/rate  expression  if  none  of  the  changes 

to  selected  results  of  the  simulations  exceed  specified  criteria. 

A  species  is  eliminated  as  a  consequence  of  all  the  reactions  involving  it  being 
eliminated. 

As  in  previous  efforts  (Kotlar  2010;  McQuaid  2013;  McQuaid  and  Chen  2014), 
parameters  of  the  simulations  that  were  compared  were  local  maxima  found  in 
mass-specific  and  volumetric  “heat  release  rate”  ( Qmass  ar|d  respectively) 

versus  time  (t)  traces,  the  times  at  which  those  maxima  occurred  (t^ass  and  t™“x, 
respectively),  and  the  temperature  at  the  final  time  step  of  the  simulation  (T/).  In 
contrast  to  the  efforts  by  Kotlar  (2010)  and  McQuaid  (2013),  however,  qmass  and 
qvoi  were  calculated  per 

RmassCt)  =  R  (m  *  Yk(t)/Wk)  (la) 

and 

Qvol(t)  Pg(t)  *  <?jnass(0>  (l^t) 

where  R  is  the  universal  gas  constant,  T  is  the  temperature,  Yk  is  the  mass  fraction 
of  the  Kh  species,  Wk  is  its  molecular  weight,  and  pg  is  the  gas  mixture’s  density. 
This  approach  was  developed  in  the  course  of  a  prior  effort  to  avoid  numerical 
instabilities  that  occasionally  arose  when  qmass  and  qvol  were  computed  on  the 
basis  of  species’  chemical  production  rates  and  their  enthalpies  (McQuaid  and  Chen 
2014).  Traced  to  poorly  formulated/”spliced”  coefficients  employed  for  calculating 
species’  enthalpies,  that  issue  has  been  resolved  and  does  not  apply  to  the  current 
AP-HTPB  mechanism.  However,  since  screening  based  on  Eqs.  la  and  lb  is  as 
effective  and  is  less  computationally  expensive  than  the  previous  method,  there  was 
no  reason  to  return  to  prior  practice. 

Following  the  same  reaction  step  order  (minus  those  that  had  been  eliminated),  the 
screening  process  was  repeated  for  a  given  set  of  acceptance  criteria  if  at  least  one 
reaction  was  permanently  eliminated  in  the  course  of  a  complete  pass  through  the 
mechanism.  When  no  more  reactions  could  be  eliminated  without  exceeding  any 
specified  tolerance,  or  10  complete  passes  at  a  specified  set  of  tolerances  had  been 
completed,  the  tolerances  were  relaxed  and  the  process  repeated. 

The  HR  model  employed  to  perform  the  simulations  used  the  CHEMKIN-II 
subroutine  library  (Kee  et  al.  1989).  Solutions  for  the  model’s  differential-algebraic 
equations  were  obtained  with  DASPK  (Li  and  Petzold  2000).  For  the  reductions 
reported  here,  the  tolerances  for  variations  in  q™“x,  q^ass ,  0^5  and  t^ass  were 
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set  at  ±1.0  %  for  the  first  pass.  The  tolerance  for  variation  in  7/ was  set  at  ±1.0  K. 
They  were  subsequently  relaxed  in  ±1.0%— ±1.0  K  increments.  The  largest 
tolerances  employed  were  ±10%  and  ±10  K.  Collectively  they  are  referred  as  the 
±10%— ±10  K  acceptance  criteria.  As  will  be  discussed,  the  bases  for  these  criteria 
are  empirical  and  not  foolproof.  Thus  postreduction  analyses  are  an  integral  aspect 
of  candidate  selection.  In  the  present  case,  they  were  needed  to  address  issues  that 
arose  because  trial  mechanisms  could  not  be  screened  (initially)  on  the  basis  of 
simulations  that  covered  the  entire  parameter  space  considered  relevant. 

Table  1  shows  the  input  parameters  for  sets  of  HR  simulations  that  were  employed 
for  the  reduction  effort.  Starting  with  set  A  (and  evolving  from  B,  to  C,  to  D  as  the 
effort  progressed),  they  were  based  on  observations  pertaining  to  AP-HTPB 
composite  propellants  and  their  deflagration. 

Table  1  Input  parameters  for  CONP  HR  simulations  employed  for  the  initial  reduction  of 
the  full  mechanism 


Set 

Start3 

Comp.b 

Sim. 

Wt%  AP 

wt% 

HTPB 

Initial 

Temp. 

(K) 

Pressure 

(psia) 

Sim. 

Length 

(S) 

A 

96 

25 

1 

100 

0 

900 

3,000 

1.0000 

A 

2 

80 

20 

900 

3,000 

5.0000 

B 

48 

>llc 

1 

100 

0 

900 

1,000 

0.2000 

2 

80 

20 

900 

1,000 

0.0100 

1 

100 

0 

900 

1,000 

0.2000 

C 

48 

29 

2 

80 

20 

900 

1,000 

0.0100 

3 

80 

20 

900 

10,000 

0.0015 

1 

100 

0 

900 

1,000 

0.2000 

D 

48 

7 

2 

80 

20 

900 

1,000 

0.0100 

/ 

3 

80 

20 

900 

10,000 

0.0015 

4 

100 

0 

900 

10,000 

0.0060 

aNumber  of  reductions  with  different  reaction  orderings  that  were  started. 

bNumber  of  different  reaction  orderings  for  which  at  least  one  pass  through  the  full  mechanism  was  completed 
and  a  workable  reduced  mechanism  was  produced. 

cSome  files  were  inadvertently  deleted,  and  the  actual  number  of  orderings  that  produced  at  least  one  reduced 
mechanism  is  no  longer  known  but  it  was  greater  than  1 1 . 


As  in  prior  efforts  (Kotlar  2010;  McQuaid  2013;  McQuaid  and  Chen  2014),  the 
starting  chemical  compositions  were  chosen  to  reflect  expectations  for  the  nascent 
products  of  the  condensed-phase  to  gas-phase  conversion  process  of  individual 
propellant  ingredients,  in  this  case  AP  and  HTPB.  For  the  nascent  product  of  AP’s 
decomposition,  an  NH3-HCIO4  complex  (labeled  NH4CIO4)  was  chosen  based  on 
prior  work  (McQuaid  and  Chen  2014).  In  that  work,  an  earlier  version  of  the  full 
mechanism  discussed  in  this  report  was  reduced  and  a  resultant  mechanism 
employed  with  the  CYCLOPS  framework  to  calculate  pure  AP’s  burning  rate  as  a 
function  of  pressure.  Finding  conflicting  evidence  as  to  the  nature  of  the  nascent 
product(s)  of  the  condensed-phase  to  gas-phase  conversion  process — with  either 
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NH4CIO4  or  NH3  +  HCIO4  being  the  primary  candidates — we  compared  burning- 
rate  estimates  produced  by  them  individually.  The  difference  was  found  to  be 
negligible.  The  decision  to  employ  NH4CIO4  instead  of  NH3  +  HCIO4  (or  some 
combination  of  the  2)  was  made  because  NH4CIO4  was  more  “condensed-phase- 
like”  and  it  was  easy  to  imagine  it  as  a  transitional  state  for  AP’s  decomposition. 
Burning-rate  estimates  that  were  in  reasonable  agreement  with  measured  values 
affirmed  the  choice. 

A  relatively  large  hydrocarbon  (C20H32)  with  the  same  ratio  of  cis/trans/vinyl 
butadiene  substructures  as  R45M  polymer  was  chosen  as  the  nascent  product  of 
HTPB’s  pyrolysis.  Labeled  R45M,  CYCLOPS-based  calculations  that  employed  it 
as  the  nascent  product  of  HTPB’s  thermally  induced  decay  were  found  to  yield  rates 
that  were  in  good  agreement  with  measured  values  (Chen  and  McQuaid  2009). 

The  choice  of  NH4C104:R45M  concentration  ratios  to  use  to  start  the  HR 
simulations  was  less  straightforward.  The  relative  weight  of  AP  in  AP-HTPB 
formulations  (without  aluminum)  is  typically  in  the  range  from  80%  to  90% 
(Cai  et  al.  2008;  Hedman  et  al.  2012,  Gross  et  al.  2013).  However,  as  mentioned  in 
the  Introduction,  the  range  of  fuel-oxidizer  ratios  in  the  gas  phase  near  the  surface 
of  a  burning  composite  propellant  is  much  broader  than  that.  Being 
composite/nonhomogeneous,  there  are  regions  near  the  surface  that  are  very  rich  in 
the  daughters  of  HTPB’s  pyrolysis  and  others  that  are  rich  in  the  daughters  of  AP’s 
decomposition  (Surzhikov  and  Krier  2003;  Gross  et  al.  2013).  As  such,  to  increase 
the  probability  that  the  method  would  produce  candidates  that  were  valid  over  the 
entire  pertinent  range,  simulations  with  starting  compositions  that  ranged  from 
1 00%  NH4CIO4  to  1 00%  R45M  would  normally  have  been  employed  for  screening. 
However,  for  the  TMM  to  be  practical,  one  complete  pass  through  the  full 
mechanism  needs  to  be  complete  in  less  than  168  h  (1  week)  of  central  processing 
unit  time,  and  that  forced  some  compromises.  The  CONP  simulations  chosen  as  the 
basis  for  the  reduction  proved  to  have  ran  times  of  2-5  min  when  the  HR  model 
was  equipped  with  the  full  mechanism,  and  the  full  mechanism  had  2,627  reactions. 
Therefore,  a  single  pass  through  the  mechanism  had  the  potential  to  take  more  than 
a  week  when  just  one  simulation  was  employed  for  the  screening  process. 
(Simulation  times  get  shorter  as  reactions  are  eliminated.  Thus  the  total  time  needed 
to  complete  a  pass  depends  on  the  number  of  reactions  that  are  eliminated,  making 
a  priori  estimates  for  the  time  it  will  take  to  complete  the  first  pass  unreliable.)  As 
such,  the  number  of  different  simulations  that  could  be  employed  as  bases  for  the 
decision  to  eliminate  a  reaction  was  severely  limited,  preventing  us  from  relying 
completely  on  the  TMM  to  screen  candidates  over  the  entire  composition  range  of 
interest.  Rather,  we  needed  to  choose  a  restricted  range  over  which  to  screen 
eliminations,  then  search  for  candidates  that  met  acceptance  criteria  outside  that 
range. 
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Given  those  considerations,  we  made  some  judgements  about  the  relative 
importance  of  fuel-oxidizer  ratios  within  the  full  range.  First,  we  concluded  that  the 
single  most  important  ratio  to  include  would  be  for  the  80-20  AP-HTPB  case. 
Having  an  fuel-oxidizer  equivalence  ratio  close  to  1 ,  this  case  would  ensure  that  the 
process  went  to  completion,  and  it  was  assumed  it  would  be  the  one  most  likely  to 
be  able  to  produce  a  reduced  mechanism  that  could  represent  the  combustion 
chemistry  over  the  entire  range  of  compositions  that  are  pertinent  to  the  process. 

Next,  observing  that  the  deflagration  of  pure  AP  is  self-sustaining  at  pressures 
above  20  atm  (300  psi)  and  that  it  is  clearly  an  important  subprocess  in  the 
deflagration  of  AP-HTPB  composite  propellants  (Gross  et  al.  2011;  Gross  et  al. 
2013),  we  considered  it  important  that  the  reduced  mechanism  be  able  to  produce 
simulations  that  mimic  those  produced  by  the  full  mechanism  for  this  composition. 
The  complete  decomposition/pyrolysis  of  100-wt%  HTPB  in  the  absence  of 
oxidizers,  on  the  other  hand,  was  assumed  to  be  far  less  important.  It  is  slow  on  the 
1-  to  10-ms  time  scale  needed  for  the  heat-up  and  conversion  of  a  propellant  into 
final  products  (Chen  and  McQuaid  2009).  Thus,  we  did  not  consider  it  necessary 
that  the  results  produced  by  a  reduced  mechanism  meet  the  same  criteria  for 
modeling  HTPB’s  pyrolysis  as  those  for  the  stoichiometric  to  oxidizer-rich 
conditions.  Based  on  these  considerations,  and  not  expecting  that  more  than  2 
simulations  could  be  employed  for  screening,  we  chose  100-wt%  NH4CIO4  as  the 
starting  composition  for  one  simulation  and  80-20-wt%  NH4CIO4-R.45M  as  the 
starting  composition  for  the  other. 

For  all  the  HR  simulations  performed  in  this  effort  the  starting  temperature  was 
chosen  to  be  900  K.  This  choice  stemmed  from  observations  made  in  the  course  of 
reducing  an  earlier  version  of  the  AP  (only)  submechanism  of  the  current  full 
AP-HTPB  mechanism  and  burning  rate  estimates  for  pure  AP  produced  by 
PREMIX/CYCLOPS-based  models  (McQuaid  and  Chen  2014).  In  those  efforts, 
pressure-dependent  starting  temperatures  ranging  from  800  to  1 ,200  K  were  chosen 
because  published  pyrolysis  laws  indicated  that  the  surface  temperature  of 
deflagrating  AP  would  be  in  this  range  and  rise  with  pressure. 

That  said,  the  correspondence  between  the  starting  conditions  for  the  HR 
simulations  and  the  conditions  employed  in  the  CYCLOPS-based  deflagration 
model  was  not  as  clear  cut  as  imagined.  Briefly,  we  have  always  assumed  that  one 
of  the  reasons  HR  simulations  are  an  effective  basis  for  producing  reduced 
mechanisms  for  use  in  modeling  combustion  chamber  dynamics  and  deflagration 
phenomena  is  that  there  is  a  fairly  direct  correspondence  between  the  temporal 
dimension  of  HR  simulations  and  the  spatial  dimension(s)  of  the  gas-phase 
combustion  of  a  burning  propellant.  However,  in  the  present  case  there  was  a  very 
rapid  equilibration  step  between  NH4CIO4  and  NH3  +  HCIO4  in  the  HR  simulation 
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that  dropped  the  temperature  more  than  100  K  within  about  10'9  s.  As  will  be 
discussed,  nothing  comparable  happens  in  the  spatial  dimension  of  the  1-D  flame 
simulations.  Nevertheless,  we  found  that  the  reduced  mechanism  derived  on  the 
basis  of  the  aforementioned  HR  simulations  yielded  reasonable  CYCLOPS-based 
estimates  for  the  linear  burning  rate  of  pure  AP.  We  also  found  that  linear  burning 
rates  were  effectively  independent  of  the  temperature  that  was  specified  for  the 
surface  (as  dictated  by  a  pyrolysis  law)  (McQuaid  and  Chen  2014).  Thus,  900  K 
was  deemed  a  reasonable  choice  to  employ  for  all  pressures. 

Not  expecting  that  more  than  2  simulations  would  be  able  to  be  employed  as  bases 
for  the  (initial)  reduction,  and  that  one  should  have  80-20-wt%  NH4CIO4-R45M  as 
the  starting  composition  and  100-wt%  NH4CIO4  for  the  other,  we  decided  to  use 
one  (intermediate)  pressure  for  both,  then  search  for  candidates  that  would  produce 
results  similar  to  the  full  mechanism  over  the  pressure  range  of  particular  interest: 
1,000-10,000  psia.  We  chose  3,000  psia  for  the  first  attempt  (set  A)  since  exponent 
breaks  generally  occur  near  that  value.  However,  that  set  did  not  produce  a 
candidate  that  was  able  to  reproduce  a  secondary  (trass’  Qmass )  datum  produced  at 
1,000  psia  (Fig.  1.)  The  pressures  specified  for  the  simulations  employed  in  sets 
B-D  addressed  this  issue. 
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Fig.  1  Comparison  of  mass-specific  heat  release  and  temperature  vs.  time  traces  for  CONP 
HR  simulations  produced  with  the  full,  Rl,  and  R2  mechanisms;  1,000  psia 
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Considering  it  important  that  simulations  approach  an  equilibrium  state,  we 
specified  temporal  durations  for  set  A  such  that  the  temperature  difference  between 
the  final  2  time  steps  of  a  simulation  was  effectively  zero.  However,  with  such 
durations  employed,  most  of  the  reduced  mechanisms  that  were  produced  generated 
T  versus  t  traces  that  for  t  near  the  primary  heat  release  event  were  by  visual 
inspection  in  poor  agreement  with  simulations  produced  by  the  full  mechanism. 
The  shorter  durations  of  the  simulations  in  sets  B-D  were  specified  to  address  this 
issue. 

Candidates  with  100  or  fewer  reactions  and  100  or  fewer  species  were  produced  by 
sets  A-C.  However,  postreduction  analyses  based  on  visual  inspection  of  full 
T  versus  t  traces  raised  concerns  about  all  of  them.  Therefore,  T  versus  t  traces 
meeting  ±1.0%-±1.0  K  acceptance  criteria  in  sets  A  and  B  were  visually  inspected 
to  identify  candidates  that  might  be  suitable  for  further  reduction  with  a  wider  set 
of  screening  criteria.  That  search  led  to  the  identification  of  a  candidate  from  set  B 
that  produced  traces  that  well-mimicked  those  produced  by  the  full  mechanism  over 
the  entire  parameter  space  of  interest.  Referred  to  hereafter  as  Rl,  it  had  203 
reactions  and  involved  123  species. 

Subsequent  to  Rl’s  identification,  a  mistake  in  the  full  mechanism  and  Rl  was 
identified  and  corrected.  This  change  had  a  negligible  impact  on  the  T  versus  t 
traces  produced  by  the  full  mechanism  and  Rl.  It  did,  however,  result  in  Rl  failing 
to  produce  2  qmass  maxima  in  the  AP-HTPB  simulation  at  1 ,000  psia.  Appearing 
as  a  shoulder  on  the  main  feature  of  the  qmass  versus  t  trace,  the  feature  associated 
with  the  second  maximum  is  clearly  observed.  Thus,  this  discrepancy  was 
considered  ignorable. 

Rl  was  reduced  based  on  the  set  of  simulations  whose  input  parameters  are  listed 
in  Table  2.  When  postreduction  analyses  indicated  a  mechanism  produced  by  its 
reduction  met  most  of  the  criteria  that  were  sought,  the  reduction  effort  was 
stopped.  That  mechanism  is  referred  to  hereafter  as  R2. 

Table  2  Input  parameters  for  CONP  HR  simulations  employed  for  the  reduction  of  Rl 


Set 

Start3 

Comp.b 

Sim. 

Wt%  AP 

wt% 

HTPB 

Initial 

temp. 

(K) 

Pressure 

(psia) 

Sim. 

length 

(S) 

1 

100 

0 

900 

500 

0.3000 

48 

48 

2 

80 

20 

900 

500 

0.0120 

J 2/ 

3 

100 

0 

900 

20,000 

0.0015 

4 

80 

20 

900 

20,000 

0.0006 

aNumber  of  reductions  with  different  reaction  orderings  that  were  started. 

bNumber  of  different  reaction  orderings  for  which  at  least  one  pass  through  the  full  mechanism  was  completed 


and  a  workable  reduced  mechanism  was  produced. 
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3.  Postreduction  Analysis 


The  criteria  employed  in  this  effort  for  deciding  whether  a  simulation  produced  by 
a  trial  mechanism  was  in  good  agreement  with  a  simulation  produced  by  the  full 
mechanism  are  to  a  large  extent  the  same  as  those  that  were  employed  by  Kotlar 
(2010)  in  the  first  published  application  of  the  TMM.  That  application  and 
subsequent  ones  have  borne  out  their  efficacy  but  they  are  not  foolproof.  Decisions 
to  eliminate  reactions  are  based  on  comparing  (for  each  simulation)  1  or  2  data 
points  from  3  different  traces  that  usually  have  in  excess  of  1 ,000  data  points  each. 
Therefore,  it  is  possible  (and  sometimes  frequently  observed)  that,  when  plotted  in 
their  entirety  and  compared  with  the  standard,  simulations  meeting  all  the 
acceptance  criteria  exhibit  traits  that  raise  concerns.  Due  to  such  “false  positives”, 
comparing  the  traces  produced  by  candidates  with  those  produced  by  the  standard 
is  a  necessary  aspect  of  candidate  validation  and  selection.  Nevertheless,  to 
determine  if  HR  simulations  produced  with  R1  and  R2  would  meet  the  acceptance 
criteria  employed  in  the  reduction  process  over  the  entire  parameter  space  of 
interest,  a  set  of  simulations  that  spanned  the  space  were  run  and  discrepancies 
identified.  The  set  included  simulations  at  pressures  from  300  to  22,000  psia 
(20-1,500  atm).  The  results  are  shown  in  Tables  3  and  4. 
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Table  3  Comparison  of  parameters  produced  by  the  full  and  reduced  mechanisms  in  CONP  HR 
simulations  in  which  the  starting  composition  was  100-wt%  NH4CIO4 


Pressure 

(psia) 

Parameter 

Units 

Full 

R1 

A 

R2 

A 

tmax  a 

Lvol/mass 

seconds 

4.123E-01 

3.970E-01 

-4% 

3.236E-01 

-21% 

•max 

Hvol 

cal/cm3-s 

1.427E+07 

1.418E+07 

-1% 

1.316E+07 

-8% 

300 

cal/g-s 

2.415E+09 

2.41 1E+09 

0% 

2.294E+09 

-5% 

?S(2) 

cal/g-s 

2.345E+09 

2.277E+09 

-3% 

2.439E+09 

+4% 

7X1  s) 

K 

2686 

2685 

-1  K 

2686 

OK 

j.max  a 

Lvol/mass 

seconds 

2.073E-01 

2.049E-01 

-1% 

2.050E-01 

-1% 

•max 

Hvol 

cal/cm3-s 

4.039E+07 

4.003E+07 

-1% 

3.693E+07 

-9% 

500 

4Z?S(  i) 

cal/g-s 

4.149E+09 

4.134E+09 

0% 

3.91 1E+09 

-6% 

«fs(2) 

cal/g-s 

3.884E+09 

3.798E+09 

-2% 

4.046E+09 

+4% 

?X1  s) 

K 

2706  K 

2705 

-1  K 

2706  K 

OK 

j.max  a 

Lvol/mass 

seconds 

7.670E-02 

7.730E-02 

+1% 

8.386E-02 

+9% 

•max 

Hvol 

cal/cm3-s 

1.557E+08 

1.540E+08 

-1% 

1.416E+08 

-9% 

1,000 

cal/g-s 

8.336E+09 

8.298E+09 

0% 

7.853E+09 

-6% 

«?s(2) 

cal/g-s 

7.490E+09 

7.453E+09 

0% 

7.863E+09 

+5% 

7>(1  s) 

K 

2730 

2730 

OK 

2730 

OK 

u.max  a 

Lvol/mass 

seconds 

1.457E-02 

1.480E-02 

+2% 

1.639E-02 

+12% 

3  000 

•max 

Hvol 

cal/cm3-s 

1.243E+09 

1.228E+09 

-1% 

1.128E+09 

-9% 

J  ^  V/  V/ 

xmax 

Hmass 

cal/g-s 

2.416E+10 

2.340E+10 

-3% 

2.274E+10 

-6% 

7X1  s) 

K 

2763 

2764 

+1  K 

2764 

+1  K 

tmax  a 

Lvol/mass 

seconds 

2.330E-03 

2.333E-03 

0% 

2.576E-03 

+11% 

1 0  000 

•max 

Hvol 

cal/cm3-s 

1.238E+10 

1.233E+10 

0% 

1.132E+10 

-9% 

xmax 

Hmass 

cal/g-s 

7.591E+10 

7.586E+10 

0% 

7.127E+10 

-6% 

7X1  s) 

K 

2794 

2796 

+2  K 

2796 

+2  K 

j-max  a 

Lvol/mass 

seconds 

9.136E-04 

9.025E-04 

-1% 

9.931E-04 

+9% 

90  000 

•max 

Hvol 

cal/cm3-s 

4.399E+10 

4.451E+10 

+  1% 

4.062E+10 

-8% 

xvj  ,vj\j\j 

xmax 

Hmass 

cal/g-s 

1.399E+1 1 

1.416E+1 1 

+  1% 

1.313E+1 1 

-6% 

7X1  s) 

K 

2810 

2812 

+2  K 

2812 

+2  K 

aTo  the  number  of  significant  digits  necessary  to  establish  whether  the  acceptance  criteria  were  met,  £™xand  t^ass  were  the 
same. 
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Table  4  Comparison  of  parameters  produced  by  the  full  and  reduced  mechanisms  in  CONP  HR 
simulations  in  which  the  starting  composition  was  80-wt%  NH4CIO4  and  20-wt%  R45M 


Pressure 

(psia) 

Parameter 

Units 

Full 

R1 

A 

R2 

A 

tmax  a 

Lvol/mass 

seconds 

4.451E-03 

4.473E-03 

0% 

4.684E-03 

+5% 

•max 

Hvol 

cal/cm3-s 

8.961E+06 

9.046E+06 

+1% 

8.238E+06 

-8% 

300 

ffis(l) 

cal/g-s 

1.286E+09 

1.301E+09 

+  1% 

1.230E+09 

-4% 

?S(2) 

cal/g-s 

1.207E+09 

1.224E+09 

+1% 

1.236E+09 

+2% 

7X1  s) 

K 

3143 

3144 

+1  K 

3162 

+19  K 

tmax  a 

Lvol/mass 

seconds 

3.796E-03 

3.829E-03 

+1% 

3.949E-03 

+4% 

•max 

Hvol 

cal/cm3-s 

2.503E+07 

2.514E+07 

0% 

2.289E+07 

-9% 

500 

cal/g-s 

2.180E+09 

2.195E+09 

+1% 

2.076E+09 

-5% 

xmax  (o) 
Hmass 

cal/g-s 

1.954E+09 

1.988E+09 

+2% 

2.033E+09 

+4% 

7X1  s) 

K 

3182 

3184 

+2  K 

3202 

+20  K 

j-max  a 

Lvol/mass 

seconds 

2.817E-03 

2.829E-03 

0% 

2.859E-03 

+1% 

•max 

Hvol 

cal/cm3-s 

9.461E+07 

9.456E+07 

0% 

8.646E+07 

-9% 

1,000 

qZassW 

cal/g-s 

4.288E+09 

4.302E+09 

0% 

4.094E+09 

-5% 

G?s(  2) 

cal/g-s 

3.666E+09 

b 

b 

7X1  S) 

K 

3228 

3230 

+2K 

3249 

+21  K 

j.max  a 

Li wl/mass 

seconds 

1.393E-03 

1.403E-03 

+  1% 

1.376E-03 

-1% 

2  non 

•max 

Hvol 

cal/cm3-s 

7.341E+08 

7.270E+08 

-1% 

6.739E+08 

-8% 

xmax 

Hmass 

cal/g-s 

1.203E+10 

1.195E+10 

-1% 

1.151  E+ 1 0 

-4% 

7X1  s) 

K 

3288 

3290 

+2  K 

3311 

+24  K 

f max  a 

Li wl/mass 

seconds 

4.885E-04 

5.021E-04 

+3% 

4.859E-04 

-1% 

1  o  non 

•max 

Hvol 

cal/cm3-s 

7.141E+09 

6.980E+09 

-2% 

6.563E+09 

-8% 

•max 

Hmass 

cal/g-s 

3.674E+10 

3.581E+10 

-3% 

3.455E+10 

-6% 

7/(1  s) 

K 

3338 

3340 

+2  K 

3364 

+26  K 

j.max  a 

Li wl/mass 

seconds 

2.653E-04 

2.752E-04 

+4% 

2.672E-04 

+1% 

?o  non 

•max 

Hvol 

cal/cm3-s 

2.501E+10 

2.437E+10 

-3% 

2.281E+10 

-9% 

xmax 

Hmass 

cal/g-s 

6.654E+10 

6.436E+10 

-3% 

6.128!pE+10 

-8% 

7/(1  s) 

K 

3358 

3361 

+3  K 

3387 

+29  K 

aTo  the  number  of  significant  digits  necessary  to  establish  whether  the  acceptance  criteria  were  met,  t^xand  t^ass 
were  the  same. 

bThe  maximum  appearing  in  the  simulation  produced  by  the  full  mechanism  appears  as  a  shoulder  that  does  not  have  a 
local  maximum  when  R1  and  R2  were  employed  to  produce  the  simulation. 


Table  3  presents  the  results  that  were  obtained  for  the  pure  100-wt%  NH4CIO4  HR 
simulations.  It  is  observed  that  only  t™°fimass  values  in  simulations  produced  by 
R2  fail  to  meet  the  ±10%— ±10  K  acceptance  criteria.  We  do  not  believe  these 
discrepancies  disqualify  R2  for  use  in  intended  applications.  Kotlar  (2010) 
uniformly  employed  ±10%  for  variations  in  q™®*,  q™assi  and  t^ass  in  his 

initial  implementation  of  the  TMM.  Reported  without  precedent  or  explanation,  this 
choice  was  presumably  based  on  his  intuition  and/or  observations  he  made  in  the 
course  of  developing  the  method.  Given  their  effectiveness/efficacy  in  that 
application  and  all  mechanism  reduction  efforts  undertaken  by  our  group  since 
then,  there  has  been  no  reason  to  critically  examine  the  impact  of  relaxing  the 
acceptance  criterion  for  specific  parameters.  (There  have,  however,  been  cases  in 
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which  alternate  criteria  have  been  employed.)  Indeed,  given  the  number  of  variables 
involved,  the  complexity  of  their  interactions,  and  the  fundamental  differences 
between  a  homogeneous  reactor  and  combustion  scenarios  involving  spatial 
dimensions,  it  would  be  quite  difficult  to  quantify. 

That  said,  observations  regarding  the  difference  in  sensitivity  to  starting  conditions 
of  t^/mass  values  in  HR  simulations  and  burning  rate  predictions  based  on  the 
PREMIX/CYCLOPS  framework  indicate  that  the  ±10%  acceptance  criterion  for 
this  parameter  is  probably  more  restrictive  than  it  needs  to  be.  For  example,  the 
difference  in  the  t^°^mass  values  produced  by  the  full  mechanism  for  the  AP 
system  at  500  and  1,000  psia  are  0.2071  and  0.0768  s,  respectively — a  270% 
difference/increase.  This  difference  is  almost  4  times  larger  than  the  70%  difference 
in  CYCLOPS-based  estimates  for  AP’s  burning  rates  at  these  pressures  that  were 
produced  with  a  reduced  mechanism  derived  from  an  earlier  version  of  the  current 
full  mechanism  (McQuaid  and  Chen  2014). 

Even  more  telling  and  relevant  is  the  sensitivity  of  t™°fimass  to  the  starting 
temperature  of  the  simulation.  For  example,  when  the  HR  model  is  equipped  with 
the  full  mechanism  and  the  starting  temperature  of  the  300  psia,  the  100-wt% 
NH4CIO4  simulation  is  reduced  1  K  (from  900  to  899  K),  t^LJmass  increases  from 
4.123  x  10"1  s  to  4.389  x  10"1  s — a  6  %  difference.  This  is  notable  because  although 
the  chemical  compositions  and  starting  temperature  of  the  simulations  that  are 
compared  for  mechanism  reduction  are  identical  at  t  =  0  s,  as  already  mentioned, 
the  model  HR  system  undergoes  an  initial  “equilibration”  step  that  significantly 
lowers  the  temperature.  In  the  100- wt%  NH4CIO4  simulation  at  300  psia,  a  drop  of 
almost  300  K  occurs  in  less  than  1  x  10~9  s.  Having  no  analog  in  the  deflagration 
model  results  discussed  in  the  next  section,  this  aspect  of  the  simulation  was  not 
considered  in  the  mechanism  screening  process.  Differences  between  results 
produced  by  the  full  and  reduced  mechanisms  for  the  process  at  this  early  time  are 
noticeable  and  undoubtedly  contribute  to  the  t^“*mass  discrepancies. 

One  could,  of  course,  attempt  to  address  this  issue  by  setting  a  tolerance  on  the 
temperature  variation  at  some  arbitrarily  selected  temporal  point  prior  to  the 
primary  heat  release  event.  However,  since  we  do  not  believe  that  a  21%  difference 
in  the  t™°fimass  values  produced  by  R2  and  the  full  mechanism  will  result  in  as 
large  a  difference  in  the  burning  rate  estimates  produced  by  them,  the  imposition  of 
such  a  criterion  would  likely  have  excluded  potentially  viable  candidates. 
Moreover,  even  if  the  maximum  difference  in  the  burning  rate  estimates  produced 
by  the  2  approached  21%,  from  the  standpoint  of  studying  the  pressure-dependent 
burning  rates  of  AP-HTPB  composite  propellants,  it  would  be  acceptable  as  long 
as  it  varied  slowly  and  consistently  with  pressure,  preventing  the  introduction  of 
artificial  trends  or  the  masking  of  real  ones.  Regardless,  since  the  final  arbiter  of 
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the  validity  of  relaxing  this  criterion  will  be  how  the  results  of  a  deflagration  model 
equipped  with  R2  compare  with  experimental  data,  this  issue  is  simply  something 
about  which  potential  users  should  be  aware. 

Table  4  shows  the  results  that  were  obtained  for  the  80-20  NH4CIO4-R45M  cases. 
It  is  observed  that  all  the  simulations  produced  by  R1  meet  the  ±10%-±10  K 
acceptance  criteria.  Except  for  Tf  values,  that  is  also  the  case  for  the  simulations 
produced  by  R2. 

The  Tf  exceptions  require  explanation.  As  mentioned  in  Section  2,  times  on  the 
order  of  1  s  were  needed  for  the  temperature  difference  between  successive  time 
points  to  approach  0,  and  that  was  the  duration  employed  for  the  postreduction 
analysis.  However,  screening  mechanisms  on  the  basis  of  the  Tf  produced  at  such 
durations  yielded  candidates  that  generated  simulations  that  did  not  well-mimic  the 
T  versus  t  trace  near  the  major  heat  release  feature.  To  address  that  issue,  a  much 
shorter  time  to  compare  “final”  temperatures  was  specified,  allowing  deviations  to 
exceed  ±10  K  at  1  s. 

Considering  this  unexpected  outcome,  we  noted  that  Kotlar  (2010)  used  (again 
without  precedent  or  explanation)  ±200  K  for  the  first  published  application  of  the 
TMM,  and  mechanisms  produced  with  it  subsequently  demonstrated  their  ability  to 
represent  the  combustion  chemistry  of  interest.  The  employment  of  ±10  K  for  the 
reduction  process  in  the  current  study  stems  from  McQuaid  (2013)  finding  that  the 
use  of  a  tighter  tolerance  for  7/variations  produced  fewer  false  positives.  Moreover, 
the  tighter  tolerance  did  not  significantly  compromise  the  TMM’s  ability  to  produce 
mechanisms  with  targeted  sizes.  Given  the  successes  achieved  with  mechanisms 
produced  by  Kotlar  using  ±200  K,  we  see  no  reason  to  be  concerned  about 
deviations  in  7/ of  up  to  ±30  K  in  simulations  in  which  the  temperature  rises  more 
than  2,000  K. 

As  the  final  step  of  this  stage  of  the  vetting  process,  T  versus  t  and  qmass  versus  t 
traces  produced  by  Rl,  R2,  and  the  full  mechanism  were  plotted  and  compared. 
Results  for  the  AP  and  80-20  AP-HTPB  systems  at  pressures  from  300  to  22,000 
psia  (20-1,500  atm)  were  examined.  Plots  produced  at  1,000  and  10,000  psia  are 
representative  and  shown  in  Figs.  1  and  2.  Overall,  good  agreement  was  observed 
between  the  HR  simulation  results  produced  by  the  3  over  the  entire  parameter 
space  considered. 
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Fig.  2  Comparison  of  mass-specific  heat  release  and  temperature  vs.  time  traces  for  CONP 
HR  simulations  produced  with  the  full  R1  and  R2  mechanisms;  10,000  psia 
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4.  Mechanism  R2 


Rate  expressions  for  the  elementary  reaction  steps  comprising  R2  are  provided  in 
Appendix  A,  which  also  provides  information  about  the  species  involved. 
Appendix  B  provides  data  for  calculating  thermodynamic  property  estimates  for 
species  in  the  mechanism.  Appendix  C  provides  a  look-up  table  for  the  molecular 
structure  of  species  in  the  mechanism. 

5.  Burning  Rate  and  Surface  Temperature  Predictions 


5.1  Computational  Methods 

As  a  final  check  of  the  validity  of  R2  for  intended  applications,  it  was  employed  in 
a  CYCLOPS  model  to  compute  burning  rates  for  pure  AP  and  80-20  AP-HTPB. 
The  CYCLOPS  framework  is  described  in  detail  elsewhere  (Miller  and  Anderson 
2000).  Briefly,  it  involves  running  laminar,  premixed  flame  simulations  and  finding 
values  for  the  mass  flux  rate  (m)  and  temperature  at  the  burner  surface/condensed- 
phase-gas-phase  interface  ( Ts  or  T°)  such  that  2  conditions  are  satisfied.  One  is  the 
conservation  of  energy  (flux  rate)  at  the  condensed-phase-gas-phase  interface. 

Aif  «-%+0  -  >r”V”)  -  4  =0’  e> 

where  Ag  is  the  thermal  conductivity  of  the  gas  mixture  at  the  surface  (in  cal/cm-s- 
K),  Fj-0  are  the  mass  fractions  of  the  (N)  nascent  products  (i)  of  the  condensed 
phase’s  decomposition,  hf0  are  their  enthalpies  (in  cal/g)  at  Ts,  Yf°°  are  the  mass 
fractions  of  condensed-phase  ingredients,  and  h^°°  are  their  enthalpies  (in  cal/g)  at 
the  starting  bulk  temperature  (T'  ). 

The  second  condition  to  be  satisfied  is  the  relationship  between  m  and  Ts.  In  all 
prior  CYCLOPS-based  studies,  the  relationship  was  established  via  a  pyrolysis  law. 

m  =  Asexp{Es/RTs ),  (4) 

where  As  and  Es  are  constants  (Miller  and  Anderson  2000).  In  the  current  study, 
however,  the  relationship  between  m  and  Ts  was  established  through  an 
evaporative/sublimative  flux  rate  expression  (Li  et  al.  1990;  Miller  1997;  Galwey 
and  Brown  1999).  For  pure  AP,  the  expression  was 

mAP  -  Z^w^nRT^kx^Aexpi-AE/RTs)  -  PX+°)  =  0,  (5) 

Where  X^°  and  X£°  are,  respectively,  the  mole  fractions  of  the  ktb  nascent  product 
on  the  condensed  phase  and  gas-phase  sides  of  the  surface,  P  is  the  (total)  pressure, 
A  is  the  pre-exponential  for  the  equilibrium  vapor  pressure,  and  AE  is  the  effective 
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activation  energy  for  dissociative  sublimation.  A  (4.42  x  1013  dynes/cm2)  and  A E 
(28,800  cal/mol)  were  specified  on  the  basis  of  the  dissociation  pressures  of  AP  that 
were  measured  by  Inami  et  al.  (1963). 

To  approximate  the  deflagration  of  an  80-20  AP-HTPB  mixture,  we  employed  the 
homogeneous  representation  of  the  gas  phase  enabled  by  PREMIX  and  required 
that  the  linear  burning  rate  of  HTPB  ( r{?TPB  =  mHTPB/  pHTPB)  equal  the  linear 
burning  rate  of  AP  ( r£p  =  rhAP/pAP).  With  this  constraint,  it  can  be  shown  that 

ri1 AP-HTPB  —  (p AP-HTPB  /  PAp)rilAP  —  0.  (6) 


The  laminar,  premixed  flame  simulations  were  performed  with  a  model  employing 
the  CHEMKIN  III  library  of  subroutines  (Kee  et  al.  1985;  Kee  et  al.  2002).  Prior 
CYCLOPS  models  developed  and  employed  by  our  group  used  CHEMKIN  II 
subroutines.  The  change  was  made  because  simulation  results  produced  by  the 
latter  for  pure  AP  proved  overly  sensitive  to  the  presence  of  reactions  involving 
carbon-containing  molecules.  (They  should  have  had  no  impact  at  all.)  The 

(dT\+® 

—  1  was  also  revised.  Prior  versions  employed  Newton 


forward  differencing  with  a  very  limited  number  of  (x,  7)  points  to  compute 


In  this  study,  numerical  integration  of 


4  ©+°  =  Jo” 


mcvx 


dx 


(?) 


was  employed  (Miller  1997).  The  volumetric  enthalpy  release  [c/(x)],  mean  gas- 
phase  heat  capacity  (cp),  and  Ag  values  were  computed  via  CHEMKIN  III 
subroutines. 


The  condensed-phase  properties  of  AP  and  HTPB  that  were  employed  for  the 
calculations  are  shown  in  Table  5.  The  data  employed  to  compute  transport 
properties  are  given  in  Appendix  D.  Solutions  for  Ts  and  mAP /ap-htpb  were  found 
on  a  grid  with  Ts  intervals  of  0.01  K  and  m  intervals  of  0.001  g/cm2-s.  This 
resolution  yielded  residuals  for  Eq.  3  whose  magnitudes  were  less  than  0.01%  of 
m  Yii  (Yi~° ht°  ~  Yrhri  and  residuals  for  Eqs.  5  or  7  whose  magnitudes  were 
less  than  1%  of  m. 
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Table  5  Condensed-phase  properties  of  AP  and  HTPB  employed  to  calculate  burning  rates 


Material 

A Hf 
(cal/g) 

p 

(g/cm3) 

AP 

-595 

1.95 

HTPB 

-30 

0.95 

80-20  AP-HTPB 

-482 

1.75 

5.2  Results 


Table  6  presents  the  (Ts,  m )  solutions  that  were  found  for  pure  AP  burning  at 
pressures  from  300  to  20,000  psia.  The  burning  rate  estimates  as  a  function  of 
pressure  are  compared  with  measured  values  and  previously  computed  estimates  in 
Fig.  3.  It  is  observed  that  the  estimates  computed  with  R2  are  slightly  higher  than 
measured  values  at  pressures  from  300  to  2,000  psia,  and  they  are  approximately 
50%  higher  than  those  estimated  with  a  reduced  mechanism  obtained  from  a  prior 
version  of  the  AP  only  subset  of  the  current  full  mechanism.  Some  of  the  latter 

(dT\+® 

—J  was  computed. 

Given  the  limitations  of  the  model  and  that  trends  are  reproduced,  we  are  satisfied 
with  the  agreement. 

Table  6  The  surface  temperature,  mass  flux  rate  and  burning  rate  of  pure  AP  computed  at 
various  pressures  by  the  R2-CYCLOPS  model 


(atm) 

Pressure 

(psi) 

Ts 

(K) 

m 

(g/cm2-s) 

n 

(cm/s) 

20 

294 

862 

0.725 

0.37 

50 

735 

908 

1.55 

0.79 

100 

1,470 

946 

2.72 

1.39 

200 

2,940 

987 

4.70 

2.41 

400 

5,880 

1,030 

8.02 

4.11 

750 

11,000 

1,071 

12.9 

6.62 

1,500 

22,000 

1,119 

21.1 

10.8 
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Fig.  3  Comparison  of  measured  and  computed  burning  rates  for  pure  AP 


The  data  in  Table  6  also  constitutes  a  basis  for  deriving  a  pyrolysis  law.  As  shown 
in  Fig.  4,  a  plot  of  ln(m)  versus  l/Ts  is  nearly  linear,  and  the  fit  to  the  data  yields 
As  =  1.5  x  106  and  Es  =  25,200  cal/mol.  As  shown  in  Table  7,  these  values  fall  in 
the  middle  of  the  range  of  values  that  have  been  proposed/employed  by  others.  We 
also  note  that  Es  is  less  than  the  value  of  A E  (28,800  cal/mol)  employed  for  Eq.  5. 
Figure  4  shows  that  they  produce  a  relationship  between  m  and  Ts  that  is  more 
similar  to  ones  proposed  by  Beckstead  et  al.  (1970)  and  Lengelle  et  al.  (1991). 
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1/TS(K) 


Fig.  4  Comparison  of  ln(m)  vs.  1/Ts  plots  corresponding  to  different  pyrolysis  laws  for 
pure  AP 


Table  7  Comparison  of  pyrolysis  laws  proposed  for  AP  self-deflagration 


Reference 

As 

(g/cm2-s) 

Es 

(cal/mol) 

Lengelle  et  al.  (1991) 

9.6  x  104 

20,000 

Beckstead  et  al.  (1970) 

3.0  x  105 

22,000 

This  work 

1.8  x  106 

25,200 

Schusser  et  al.  (2002) 

1.6  x  108 

32,000 

Figure  5  presents  a  comparison  of  the  R2-CYCLOPS-based  burning  rate  estimates 
for  the  homogeneous  80-20  AP-HTPB  mixture  postulated  in  Section  2  with 
measured  values  for  an  80-20  AP-HTPB  mixture  formulated  with  5p  AP  particles 
(Lengelle  et  al.  2002).  (The  mixture  is  probably  as  close  to  homogeneous  as  one 
can  hope  to  attain.)  While  the  agreement  between  traces  is  not  compelling,  the  rates 
at  a  given  pressure  are  within  a  factor  of  2  of  one  another  over  the  entire  range 
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considered.  Given  the  limitations  of  the  (1-D,  laminar)  model — for  example,  for 
HTPB  to  regress  (in  nature)  at  the  same  linear  rate  as  AP,  the  temperature  of  its 
surface  must  be  much  higher  than  that  of  AP’s  surface — the  comparison  is 
considered  evidence  that  R2  will  provide  a  reasonable  representation  of  the 
gas-phase  reaction  kinetics  of  AP-HTPB  systems  in  more  detailed  CFD  models. 


Fig.  5  Comparison  of  the  pressure-dependent  burning  rates  computed  for  a  homogeneous 
80-20  AP-HTPB  mixture  and  measured  values  for  an  80-20  AP-HTPB  mixture  formulated 
with  5p  AP  particles 
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6.  Additional  Comments 


Given  the  limited  number  of  simulations  that  could  be  employed  for  the  mechanism 
reduction  process,  mechanisms  were  not  screened  for  their  ability  to  mimic  the 
oxidizer-free  pyrolysis  of  HTPB.  Based  on  the  expectation  that  reactions  involving 
oxidizers  would  dominate  the  latter  stages  of  HTPB’ s  decomposition,  this  decision 
left  open  the  possibility  that  the  reduction  process  would  produce  candidates  that 
would  not  have  a  “hydrocarbon  cracking  submechanism”  for  representing  the 
earliest  stages  of  HTPB’ s  decomposition.  However,  inspection  of  R2  indicates  that 
it  contains  a  set  of  reactions  that  could  be  considered  such  a  submechanism.  It 
includes  a  reaction  in  which  R45M  is  “cracked”  by  collision  with  any  species  in  the 
mechanism,  and  the  daughters  of  that  reaction  have  the  potential  to  decompose  via 
a  sequence  of  steps  that  do  not  involve  an  oxidizer.  Thus,  the  mechanism  appears 
to  have  the  potential  to  reasonably  represent  HTPB  pyrolysis  in  CFD  models  of 
AP-HTPB  composite  propellant  deflagration. 

7.  Summary  and  Conclusions 

A  (full)  detailed,  gas-phase,  finite-rate  chemical  kinetics  mechanism  for 
representing  the  combustion  chemistry  associated  AP-HTPB  composite  propellant 
deflagration  was  reduced  to  skeletal  sets  via  the  TMM.  The  full  mechanism  had 
2,627  elementary  reaction  steps  and  involved  637  species.  The  primary  objective 
of  the  effort  was  to  produce  a  mechanism  with  <100  reactions  and  involving  <100 
species  that  could  be  employed  as  a  submodel  in  the  ARL  CFD  laminate 
framework.  Results  produced  by  trial  mechanism-homogeneous  reactor  model 
combinations  were  the  basis  for  eliminating  reactions.  Postreduction  analyses  were 
conducted  to  identify  candidates  with  the  best  potential  to  mimic  the  full 
mechanism  in  a  parameter  space  that  is  relevant  to  AP-HTPB  composites  burning 
at  pressures  from  300  to  20,000  psia.  The  candidate  considered  to  be  the  best 
compromise  between  size  and  mimicking  ability  for  this  application  had  100 
reactions  and  involved  81  species.  Coupled  with  the  CYCLOPS  framework,  it 
produced  estimates  for  the  pressure-dependent  burning  rates  of  AP  and  a 
(hypothetical)  homogeneous  AP-HTPB  mixture  that  were  in  reasonable  agreement 
with  measured  burning  rate  data  published  in  the  open  literature.  All  data 
comprising  the  mechanism  are  provided. 
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Appendix  A.  Mechanism  R2's  Reactions,  Rate  Expressions,  and 
Species  Data 
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Table  A-l  provides  the  elementary  reaction  steps  comprising  mechanism  R2  and 
parameters  for  computing  the  rates  ( k )  of  the  steps  as  a  function  of  temperature  (7) 
and  pressure  ( P ).  Except  for  2  reactions  types  that  employ  nonstandard  methods  for 
computing  the  rates  of  fall-off  reactions,  the  reaction  descriptions  conform  to 
standard  Chemkin  II  notation  and  usage  (Kee  et  al.  1989;  references  at  end  of 
Appendix  A).  The  nonstandard  methods  for  computing  the  rates  for  fall-off 
reactions  are  designated  by  the  strings  “T&H”  (Tsang  and  Herron  1991)  and  “I nrC\ 
where  nn  is  a  number.  For  rates  computed  with  T&H  coefficients,  “Fce„?”  (Kee  et 
al.  1989)  is  calculated  per 

Fcent  =  cto  +  aiT  +  a2  T2, 

where  the  polynomial’s  coefficients  (at)  are  given  in  order  after  the  T&H  string. 
The  second  and/or  third  terms  are  zero  if  the  coefficients  are  not  provided. 

For  Inn  fall-off  reactions,  reaction  rates  are  calculated  by  logarithmically 
interpolating  between  rates  at  pressures  that  bound  the  one  of  interest  (Cantera 
2015).  That  is,  given  rate  expressions, 

k-t  (T)  =  A1Tblexp(E1/RT) 
and 

k2(J)  =  A2Tb*exp(E2/RT ) 

at  pressures  Pi  and  P2,  respectively,  the  rate  at  an  intermediate  pressure  (Pi  <  P  < 
P2)  is  computed  per 

log k(T,P)  =  logger)  +  (log k2(T)  -  logfci (7’))^Z^- 

The  coefficients  following  Inn  in  the  tables  correspond  to  At,  bi,  Ej,  and  Pi  (in 
atmospheres  [atm]),  respectively.  We  note  that  the  numbers  on  the  line  that  defines 
such  reactions  are  dummy  variables;  they  serve  no  function  in  computing  rates.  We 
also  note  that  rate  expressions  applicable  for  pressures  as  low  as  0.001  atm  were 
developed.  However,  because  R2  has  been  validated  only  for  pressures  as  low  as 
20  atm,  to  reduce  the  computational  expense  of  calculating  the  rates  for  such 
reactions,  expressions  for  pressures  less  than  10  atm  were  omitted.  That,  in  turn, 
resulted  in  certain  sets  in  the  table  that  do  not  exhibit  fall-off  behavior  in  the 
pressure  range  that  is  included.  Such  sets  could  be  written  more  economically  as  a 
simple  pressure-dependent  form  (“+M”)  with  a  single  set  of  coefficients  rather  than 
in  fall-off  form  [“(+M)”].  However,  to  prevent  such  expressions  from  being 
inadvertently  adopted  in  this  form  for  use  in  lower  pressure  applications,  the 
fall-off  form  is  specified. 

The  table  also  provides  information  about  the  species  involved,  including  their 
elemental  composition,  molecular  weight,  and  the  temperature  range  over  which 
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estimates  for  their  thermodynamic  properties  have  been  computed.  Data  for 
computing  the  species’  thermodynamic  properties  as  a  function  of  temperature  are 
provided  in  Appendix  B.  Also,  because  Chemkin  II  limits  the  number  of  characters 
that  can  be  employed  to  label  a  species  to  <16  characters,  for  many  species  it  was 
not  possible  to  assign  a  label  or  name  that,  within  the  context  of  the  mechanism 
alone,  could  be  reliably  translated  into  a  molecular  structure.  Appendix  C  provides 
a  look-up  table  that  performs  that  function. 
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Table  A-l  Species  and  reaction  rate  expressions  comprising  mechanism  R2 


ELEMENTS 

ATOMIC 

CONSIDERED 

WEIGHT 

1.  C 

12 . 0112 

2.  H 

1.00797 

3.  0 

15.9994 

4  .  N 

14.0067 

5.  CL 

35.4530 

c 

P  H 
H  A 
A  R 

SPECIES  S  G  MOLECULAR  TEMPERATURE  ELEMENT  COUNT 


CONSIDERED 

E 

E 

WEIGHT 

LOW 

HIGH 

C 

H 

O 

N 

CL 

1. 

CHO 

G 

0 

29.01852 

300 

5000 

1 

1 

1 

0 

0 

2  . 

CH20 

G 

0 

30.02649 

300 

3500 

1 

2 

1 

0 

0 

3. 

CH3 

G 

0 

15.03506 

300 

5000 

1 

3 

0 

0 

0 

4  . 

CO 

G 

0 

28.01055 

300 

5000 

1 

0 

1 

0 

0 

5. 

C02 

G 

0 

44.00995 

300 

5000 

1 

0 

2 

0 

0 

6 . 

C2H2 

G 

0 

26.03824 

300 

5000 

2 

2 

0 

0 

0 

7  . 

C2H3 

G 

0 

27.04621 

300 

5000 

2 

3 

0 

0 

0 

8. 

C*CCL*0 

G 

0 

90.50976 

300 

5000 

3 

3 

1 

0 

1 

9. 

VCD  JO 

G 

0 

55.05676 

300 

5000 

3 

3 

1 

0 

0 

10. 

VC  DO 

G 

0 

56.06473 

300 

5000 

3 

4 

1 

0 

0 

11. 

CTCV 

G 

0 

52.07648 

300 

5000 

4 

4 

0 

0 

0 

12  . 

CDCCDCJ 

G 

0 

53.08445 

300 

5000 

4 

5 

0 

0 

0 

13. 

CDCCDC 

G 

0 

54.09242 

300 

5000 

4 

6 

0 

0 

0 

14  . 

CCDOV 

G 

0 

70.09182 

300 

5000 

4 

6 

1 

0 

0 

15. 

CDCCJC 

G 

0 

55.10039 

300 

5000 

4 

7 

0 

0 

0 

16. 

VCOJC 

G 

0 

71.09979 

300 

5000 

4 

7 

1 

0 

0 

17  . 

CJW 

G 

0 

67.11154 

300 

5000 

5 

7 

0 

0 

0 

18. 

WCOJ 

G 

0 

83.11094 

300 

5000 

5 

7 

1 

0 

0 

19. 

CW 

G 

0 

68.11951 

300 

5000 

5 

8 

0 

0 

0 

20. 

CJCVC 

G 

0 

69.12748 

300 

5000 

5 

9 

0 

0 

0 

21. 

CVCCVCL 

G 

0 

144 . 64581 

300 

5000 

8 

13 

0 

0 

1 

22  . 

VY6DE13J 

G 

0 

105.16093 

300 

5000 

8 

9 

0 

0 

0 

23. 

VWJV 

G 

0 

105.16093 

300 

5000 

8 

9 

0 

0 

0 

24  . 

WVIV 

G 

0 

106.16890 

300 

5000 

8 

10 

0 

0 

0 

25. 

VY6DE13 

G 

0 

106.16890 

300 

5000 

8 

10 

0 

0 

0 

26. 

CVCZCVCJ 

G 

0 

109.19281 

300 

5000 

8 

13 

0 

0 

0 

27  . 

CVCCCOJV 

G 

0 

125.19221 

300 

5000 

8 

13 

1 

0 

0 

28. 

CJCVCCVC 

G 

0 

123.21990 

300 

5000 

9 

15 

0 

0 

0 

29. 

CVCCVCCOJ 

G 

0 

139.21930 

300 

5000 

9 

15 

1 

0 

0 

30. 

BZBCCLV 

G 

0 

198.73823 

300 

5000 

12 

19 

0 

0 

1 

31. 

BZBBJ 

G 

0 

163.28523 

300 

5000 

12 

19 

0 

0 

0 

32. 

BZBCCOJV 

G 

0 

179.28463 

300 

5000 

12 

19 

1 

0 

0 

33. 

BBBCJV 

G 

0 

203.35056 

300 

5000 

15 

23 

0 

0 

0 

34  . 

BBBCOJV 

G 

0 

219.34996 

300 

5000 

15 

23 

1 

0 

0 

35. 

BBBJVIV 

G 

0 

215.36172 

300 

5000 

16 

23 

0 

0 

0 

36. 

BBBVIV 

G 

0 

216.36969 

300 

5000 

16 

24 

0 

0 

0 

37  . 

R45MLSJ 

G 

0 

307 . 93105 

300 

5000 

20 

32 

0 

0 

1 

38. 

R45MLPJ 

G 

0 

307 . 93105 

300 

5000 

20 

32 

0 

0 

1 

39. 

R45M9J 

G 

0 

271.47008 

300 

5000 

20 

31 

0 

0 

0 

40. 

R45M13J 

G 

0 

271.47008 

300 

5000 

20 

31 

0 

0 

0 

41. 

R4  5M 

G 

0 

272.47805 

300 

5000 

20 

32 

0 

0 

0 

42  . 

CL 

G 

0 

35.45300 

300 

5000 

0 

0 

0 

0 

1 
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43. 

HCL 

G 

0 

36.46097 

300 

5000 

0 

1 

0 

0 

1 

44  . 

CLOH 

G 

0 

52.46037 

300 

5000 

0 

1 

1 

0 

1 

45. 

HOCLO 

G 

0 

68.45977 

300 

4500 

0 

1 

2 

0 

1 

46. 

HOOCL 

G 

0 

68.45977 

300 

5000 

0 

1 

2 

0 

1 

47  . 

H0CL02 

G 

0 

84 . 45917 

300 

4500 

0 

1 

3 

0 

1 

48. 

HCL04 

G 

0 

100.45857 

300 

4500 

0 

1 

4 

0 

1 

49. 

NH4CL04 

G 

0 

117 .48918 

300 

4500 

0 

4 

4 

1 

1 

50. 

AP_NH3 

G 

0 

134.51979 

300 

5000 

0 

7 

4 

2 

1 

51. 

NOCL 

G 

0 

65.45910 

300 

5000 

0 

0 

1 

1 

1 

52  . 

0N0CL03 

G 

0 

129.45670 

300 

4500 

0 

0 

5 

1 

1 

53. 

CLO 

G 

0 

51 . 45240 

300 

4500 

0 

0 

1 

0 

1 

54  . 

CL02 

G 

0 

67.45180 

300 

5000 

0 

0 

2 

0 

1 

55. 

CLOO 

G 

0 

67.45180 

300 

5000 

0 

0 

2 

0 

1 

56. 

CL03 

G 

0 

83.45120 

300 

5000 

0 

0 

3 

0 

1 

57  . 

CL2 

G 

0 

70.90600 

300 

5000 

0 

0 

0 

0 

2 

58. 

H 

G 

0 

1 . 00797 

300 

5000 

0 

1 

0 

0 

0 

59. 

HNO 

G 

0 

31.01407 

300 

5000 

0 

1 

1 

1 

0 

60. 

HONO 

G 

0 

47 . 01347 

300 

5000 

0 

1 

2 

1 

0 

61. 

NNH 

G 

0 

29.02137 

250 

4000 

0 

1 

0 

2 

0 

62. 

HNJN02 

G 

0 

61 . 02017 

300 

5000 

0 

1 

2 

2 

0 

63. 

HN0JN02 

G 

0 

77.01957 

300 

5000 

0 

1 

3 

2 

0 

64  . 

OH 

G 

0 

17 . 00737 

300 

5000 

0 

1 

1 

0 

0 

65. 

H02 

G 

0 

33.00677 

300 

5000 

0 

1 

2 

0 

0 

6  6 . 

H2 

G 

0 

2 . 01594 

300 

5000 

0 

2 

0 

0 

0 

67  . 

NH2 

G 

0 

16.02264 

300 

5000 

0 

2 

0 

1 

0 

68. 

NH20 

G 

0 

32.02204 

300 

5000 

0 

2 

1 

1 

0 

69. 

N2H2 

G 

0 

30.02934 

300 

5000 

0 

2 

0 

2 

0 

70. 

H2NN02 

G 

0 

62 . 02814 

300 

4500 

0 

2 

2 

2 

0 

71. 

H20 

G 

0 

18.01534 

300 

5000 

0 

2 

1 

0 

0 

72  . 

NH3 

G 

0 

17 . 03061 

300 

5000 

0 

3 

0 

1 

0 

73. 

N2H3 

G 

0 

31.03731 

300 

5000 

0 

3 

0 

2 

0 

74  . 

OJHNH3 

G 

0 

34.03798 

300 

5000 

0 

4 

1 

1 

0 

75. 

N2H4 

G 

0 

32.04528 

300 

5000 

0 

4 

0 

2 

0 

76. 

DINH3 

G 

0 

34 .06122 

300 

5000 

0 

6 

0 

2 

0 

77  . 

NO 

G 

0 

30.00610 

300 

5000 

0 

0 

1 

1 

0 

78. 

N02 

G 

0 

46.00550 

300 

5000 

0 

0 

2 

1 

0 

79. 

N2 

G 

0 

28.01340 

300 

5000 

0 

0 

0 

2 

0 

80. 

0 

G 

0 

15.99940 

300 

5000 

0 

0 

1 

0 

0 

81. 

02 

G 

0 

31.99880 

300 

5000 

0 

0 

2 

0 

0 

(k  =  A  T**b  exp(-E/RT))  ref. 

REACTIONS  CONSIDERED  A  b  E 


1. 

R45M (+M) =BZBB J+CVCZCVC J (+M) 

105:  1 . 9800E+91  -2.3800E+01 

124:  1 . 9800E+91  -2.3800E+01 

7 . 9841E+04 

7 . 9841E+04 

1.000E+00 

1 . 0000E+01 

2 . 0000E+03 

0.00 

o 

o 

a 

2  . 

BZBBJ (+M) =CVCZCVC J+CDCCDC (+M) 
105:  1 . 0300E+37  -7.8000E+00 

124:  1 . 0300E+37  -7.8000E+00 

3. 6423E+04 

3 . 6423E+04 

1 . 000E  +  00 

1.0000E+01 

2 . 0000E  +  03 

0.00 

0.0 

b 

3. 

VWIV=VY6DE13 

2.77  OE  +  13 

0.00 

45570.0 

c 

4  . 

VY6DE13+H=VY6DE13J+H2 

4 . 800E  +  08 

1.50 

0.0 

d 

5. 

VY  6  DE 1 3  J=VW  JV 

1.429E+12 

0.63 

61604.0 

g 

6 . 

CTCV+CDCCDC  J=VWJV 

3 . 160E  +  11 

0.00 

6000.0 

e 

7  . 

CTCV=C2H2+C2H2 

1.290E+15 

0.00 

82470.0 

f 

8. 

C2H3+C2H2=CDCCDCJ 

1.260E+10 

1.90 

2110 . 0 

h 

9. 

R4  5M+CL=R4  5M9 J+HCL 

2 . 950E+14 

0.00 

180.0 

j 

10. 

R4 5M+CLO=R4 5M1 3  J+CLOH 

3 . 443E  +  01 

3.37 

1728 . 4 

k 

11. 

R4  5M+CL=R4  5MLP J 

1 . 7  80E+10 

1.31 

-1030.0 

1 

12  . 

R4  5M+CL=R4  5MLS J 

1 . 630E+14 

0.00 

0.0 

m 

13. 

R4  5MLS J=BZBCCLV+CVCZCVC J 

1.07  6E  +  12 

0.44 

19407.0 

n 

14  . 

R4  5MLP J=BZBB J+CVCCVCL 

6 . 088E  +  11 

0.51 

18815.0 

o 

15. 

CVCZCVC J+H02=CVCCCO JV+OH 

9. 640E+12 

0.00 

0.0 

i 
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16.  WCO J=VCDO+C2H3  2.306E+13  0.30  19668.0 

17.  CDCCJC+CL02 (+M) =CCDOV+CLOH (+M)  1.000E+00  0.00  0.0 

105:  3 . 7800E+05  2.0600E+00  -1.6540E+03  1.0000E+01 

106:  3 . 7 600E  +  05  2.0600E  +  00  -1.6550E  +  03  2.5000E  +  01 

107:  3 . 8 600E  +  05  2.0600E  +  00  -1.6490E+03  5.0000E  +  01 

108:  4 . 2700E+05  2.0500E+00  -1.6260E+03  7.5000E+01 

109:  5 . 0800E+05  2.0200E+00  -1.5840E+03  1.0000E+02 

110:  8 . 2300E+05  1.9600E+00  -1.4670E+03  1.5000E+02 

Ill:  1 . 4200E+06  1.9000E+00  -1.3300E+03  2.0000E+02 

112:  2 . 4 600E  +  0 6  1.8300E  +  00  -1.1910E  +  03  2.5000E  +  02 

113:  4 . 1300E  +  0 6  1.7600E  +  00  -1.0550E  +  03  3.0000E  +  02 

114:  6 . 6900E  +  0 6  1.7000E  +  00  -9.2500E  +  02  3.5000E  +  02 

115:  1 . 0400E+07  1.6500E+00  -8.0400E+02  4.0000E+02 

116:  1 . 5600E+07  1.6000E+00  -6.9000E+02  4.5000E+02 

117:  2 . 2 600E  +  07  1.5600E  +  00  -5.8400E  +  02  5.0000E  +  02 

118:  4 . 3000E+07  1.4800E+00  -3.9100E+02  6.0000E+02 

119:  7 . 3800E+07  1.4100E+00  -2.2300E+02  7.0000E+02 

120:  1 . 1 600E+08  1.3600E+00  -7.4000E+01  8.0000E+02 

121:  1 . 7000E+08  1.3100E+00  5.8000E+01  9.0000E+02 

122:  2 . 3600E+08  1.2700E+00  1.7600E+02  1.0000E+03 

123:  6 . 6200E+08  1.1600E+00  6.1800E+02  1.5000E+03 

124:  1 . 0500E  +  0  9  1.1100E  +  00  9.0600E  +  02  2.0000E  +  03 

18.  CDCCDC=C2H3+C2H3  4.027E+19  -1.00  98150.0  q 

19.  R4  5M9  J=CW+BBBC  JV  4.480E  +  10  0.56  21910.0  r 

20.  CDCC JC+BBBVIV=R4  5M13 J  2.410E  +  10  2.48  6130.0  s 

21.  BBBCJV+N02 (+M) =BBBCOJV+NO (+M)  1.000E+00  0.00  0.0  t 

105:  9 . 6100E+20  -2.2200E+00  1.1470E+04  1.0000E+01 

106:  3 . 2500E+1 9  -1.7400E+00  1.1702E+04  2.5000E+01 

107:  8 . 5000E+17  -1.2400E+00  1.1695E+04  5.0000E+01 

108:  6 . 5400E+16  -9.0000E-01  1.1616E+04  7.5000E+01 

109:  8 . 7 600E  +  15  -6.3000E-01  1.1526E  +  04  1.0000E  +  02 

110:  3 . 9500E+14  -2.1000E-01  1.1351E+04  1.5000E+02 

Ill:  3 . 6300E+13  1.0000E-01  1.1192E+04  2.0000E+02 

112:  5 . 1 600E+12  3.6000E-01  1.1050E+04  2.5000E+02 

113:  9 . 7  900E  +  11  5.8000E-01  1.0922E  +  04  3.0000E  +  02 

114:  2 . 3000E+11  7.7000E-01  1.0804E+04  3.5000E+02 

115:  6 . 3500E+10  9.4000E-01  1.0697E+04  4.0000E+02 

116:  1 . 9900E+10  1.0900E+00  1.0597E+04  4.5000E+02 

117:  6 . 9100E  +  0  9  1.2300E  +  00  1.0504E  +  04  5.0000E  +  02 

118:  1 . 07  00E  +  0  9  1.4800E  +  00  1.0334E  +  04  6.0000E  +  02 

119:  2 . 1100E+08  1.6900E+00  1.0183E+04  7.0000E+02 

120:  5 . 0300E+07  1.8800E+00  1.0046E+04  8.0000E+02 

121:  1 . 3900E+07  2.0400E+00  9.9210E+03  9.0000E+02 

122:  4 . 3400E+06  2.2000E+00  9.8060E+03  1.0000E+03 

123:  4 . 2800E+04  2.8000E+00  9.3340E+03  1.5000E+03 

124:  1.4400E+03  3.2400E+00  8.9720E+03  2.0000E+03 

22.  BBBCO JV=BZBB J+VCDO  4.287E+11  0.52  4984.3  u 

23.  BZBBJ+N02 (+M) =BZBCCOJV+NO (+M)  1.000E+00  0.00  0.0  t 

105:  1 . 6600E+21  -2.2600E+00  1.2103E+04  1.0000E+01 

106:  1 . 9200E+19  -1.6500E+00  1.2158E+04  2.5000E+01 

107:  2 . 1 600E+17  -1.0400E+00  1.2009E+04  5.0000E+01 

108:  1 . 0200E+16  -6.3000E-01  1.1845E+04  7.5000E+01 

109:  9 . 7000E+14  -3.2000E-01  1.1696E+04  1.0000E+02 

110:  2 . 7400E+13  1.5000E-01  1.1438E+04  1.5000E+02 

Ill:  1 . 8400E+12  5.1000E-01  1.1223E+04  2.0000E+02 

112:  2 . 0500E+11  8.0000E-01  1.1038E+04  2.5000E+02 

113:  3 . 2300E+10  1.0400E+00  1.0875E+04  3.0000E+02 

114:  6 . 5000E  +  0  9  1.2500E  +  00  1.0729E  +  04  3.5000E  +  02 

115:  1 . 57  00E  +  0  9  1.4400E  +  00  1.0596E  +  04  4.0000E  +  02 

116:  4 . 4100E+08  1.6000E+00  1.0475E+04  4.5000E+02 

117:  1 . 3900E+08  1.7500E+00  1.0364E+04  5.0000E+02 

118:  1 . 8100E+07  2.0200E+00  1.0163E+04  6.0000E+02 

119:  3 . 1400E+06  2.2500E+00  9.9850E+03  7.0000E+02 
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120:  6 . 6800E+05  2.4500E+00  9.8260E+03  8.0000E+02 
121:  1 . 6800E+05  2.6300E+00  9.6820E+03  9.0000E+02 
122:  4 . 8100E+04  2.7900E+00  9.5500E+03  1.0000E+03 
123:  3 . 4800E+02  3.4300E+00  9.0150E+03  1.5000E+03 
124:  9 . 5300E+00  3.8900E+00  8.6110E+03  2.0000E+03 

24.  BZBCCO JV=C JCVCCVC+VCDO  3.451E+13  0.14  10168.0  v 

25.  CJCVCCVC+N02 (+M) =CVCCVCCO J+NO (+M)  1.000E+00  0.00  0.0  w 

105:  2 . 6100E+13  -3.0000E-02  3.5000E+01  1.0000E+01 
106:  2 . 8700E+13  -5.0000E-02  5.9000E+01  2.5000E+01 
107:  4 . 2400E+13  -9.0000E-02  1.6200E+02  5.0000E+01 
108:  6 . 8300E+13  -1.5000E-01  2.9400E+02  7.5000E+01 
109:  1 . 0800E+14  -2.1000E-01  4.2500E+02  1.0000E+02 
110:  2 . 3800E+14  -3.1000E-01  6.6400E+02  1.5000E+02 
Ill:  4 . 4300E+14  -3.8000E-01  8.6600E+02  2.0000E+02 
112:  7 . 2300E+14  -4.4000E-01  1.0370E+03  2.5000E+02 
113:  1 . 0700E+15  -4.9000E-01  1.1830E+03  3.0000E+02 
114:  1 . 4 600E  +  15  -5.2000E-01  1.3090E  +  03  3.5000E  +  02 
115:  1 . 8800E+15  -5.5000E-01  1.4190E+03  4.0000E+02 
116:  2 . 3000E+15  -5.8000E-01  1.5170E+03  4.5000E+02 
117:  2 . 7200E+15  -6.0000E-01  1.6030E+03  5.0000E+02 
118:  3 . 4 900E  +  15  -6.2000E-01  1.7500E  +  03  6.0000E  +  02 
119:  4 . 0 900E  +  15  -6.4000E-01  1.8710E  +  03  7.0000E  +  02 
120:  4 . 5200E+15  -6.5000E-01  1.9720E+03  8.0000E+02 
121:  4 . 7 600E  +  15  -6.5000E-01  2.0570E  +  03  9.0000E  +  02 
122:  4 . 8500E+15  -6.5000E-01  2.1300E+03  1.0000E+03 
123:  3 . 9500E+15  -6.1000E-01  2.3780E+03  1.5000E+03 
124:  2 . 5500E+15  -5.4000E-01  2.5150E+03  2.0000E+03 


26. 

CVCCVCCO J=CVCZCVC J+CH20 

1 . 502E+11 

0.54 

9059 

.  9 

X 

27  . 

C  VC  C  C  O  JV= C  JC  VC + VC  D  0 

3 . 451E  +  13 

0.14 

10168 

.0 

V 

28. 

BBBVIV+OH=BBB JVIV+H20 

2 . 400E  +  06 

2.00 

0 

.0 

d 

29. 

BBB  JVIV=CVCZCVC  J+WVIV 

4.019E+09 

0.86 

27888 

.0 

Y 

30. 

CO+C2H3=VCD JO 

1 . 510E+11 

0.00 

4810 

.0 

z 

31. 

C  JCVC+0=CW+0H 

3.170E+12 

0.03 

-390 

.0 

aa 

32. 

CW+OH=C  JW+H20 

3. 600E  +  06 

2.00 

0 

.0 

d 

33. 

CJW+N02  (+M)  =WCO  J+NO  (+M) 

1 . 000E+00 

0.00 

0 

.0 

w 

105:  2 . 6100E+13  -3.0000E-02  3.5000E+01  1.0000E+01 
106:  2 . 8700E+13  -5.0000E-02  5.9000E+01  2.5000E+01 
107:  4 . 2400E+13  -9.0000E-02  1.6200E+02  5.0000E+01 
108:  6 . 8300E+13  -1.5000E-01  2.9400E+02  7.5000E+01 
109:  1 . 0800E+14  -2.1000E-01  4.2500E+02  1.0000E+02 
110:  2 . 3800E+14  -3.1000E-01  6.6400E+02  1.5000E+02 
Ill:  4 . 4300E+14  -3.8000E-01  8.6600E+02  2.0000E+02 
112:  7 . 2300E+14  -4.4000E-01  1.0370E+03  2.5000E+02 
113:  1 . 0700E+15  -4.9000E-01  1.1830E+03  3.0000E+02 
114:  1 . 4 600E  +  15  -5.2000E-01  1.3090E  +  03  3.5000E  +  02 
115:  1 . 8800E+15  -5.5000E-01  1.4190E+03  4.0000E+02 
116:  2 . 3000E+15  -5.8000E-01  1.5170E+03  4.5000E+02 
117:  2 . 7200E+15  -6.0000E-01  1.6030E+03  5.0000E+02 
118:  3 . 4 900E  +  15  -6.2000E-01  1.7500E  +  03  6.0000E  +  02 
119:  4 . 0 900E  +  15  -6.4000E-01  1.8710E  +  03  7.0000E  +  02 
120:  4 . 5200E+15  -6.5000E-01  1.9720E+03  8.0000E+02 
121:  4 . 7 600E  +  15  -6.5000E-01  2.0570E  +  03  9.0000E  +  02 
122:  4 . 8500E+15  -6.5000E-01  2.1300E+03  1.0000E+03 
123:  3 . 9500E+15  -6.1000E-01  2.3780E+03  1.5000E+03 
124:  2 . 5500E+15  -5.4000E-01  2.5150E+03  2.0000E+03 


34  . 

VCDO+OH=VCD J0+H20 

1 . 000E+13 

0.00 

o 

o 

bb 

35. 

VCO JC=VCDO+CH3 

1 . 453E  +  12 

0.52 

12665.0 

cc 

36. 

VCO JC=H+CCDOV 

1.164E+10 

1.09 

12301.0 

dd 

37  . 

H+NO (+M) =HNO (+M) 

1.520E+15 

5 - 1 

o 

i 

o 

o 

ee 

Low  pressure  limit:  0.40000E+21  -0.17500E+01  0.00000E+00 


N20 

Enhanced 

by 

5 . 000E+00 

H20 

Enhanced 

by 

5 . 000E+00 

N2 

Enhanced 

by 

1 . 000E+00 
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38. 

NO+OH (+M) =HONO (+M) 

1. 

988E  +  12 

-0.05 

-721.0 

ee 

Low  pressure  limit:  0.50800E+24 

-0 . 25100E  +  01 

-0.67  600E  +  02 

T&H  VALUE:  0.62000E+00 

N20  Enhanced  by 

5 . 000E  +  00 

H20  Enhanced  by 

8 . 300E  +  00 

N2 

Enhanced  by 

1 . 000E  +  00 

39. 

HN0+N02=H0N0+N0 

4  . 

420E  +  04 

2 . 64 

4042 . 0 

ee 

40. 

HNO+O 

=OH+NO 

3. 

610E  +  13 

0.00 

0.0 

ee 

41 . 

OH+H2= 

=H20+H 

2  . 

160E  +  08 

1.50 

3430.0 

ee 

42  . 

NH2+NH2 (+M) =N2H4 (+M) 

1. 

000E  +  00 

0.00 

0.0 

ee 

105 

3 . 3100E  +  27  -4 . 5500E  +  00 

3. 9380E+03 

1 

.  0000E  +  01 

106 

4 . 7100E  +  25  -3 . 9100E  +  00 

3 . 54  90E  +  03 

2 

.  5000E  +  01 

107 

1 . 2700E  +  24  -3 . 3900E  +  00 

3 . 1 64 0E  +  03 

5 

.  0000E  +  01 

108 

1 . 4200E  +  23  -3 . 0800E  +  00 

2 . 9150E+03 

7 

.  5000E  +  01 

109 

2.9900E+22  -2.8600E+00 

2 . 7  32 0E  +  03 

1 

.  0000E  +  02 

110 

3 . 3600E+21  -2 . 5500E+00 

2 . 4680E  +  03 

1 

.  5000E  +  02 

Ill 

7 . 3300E  +  20  -2 . 3400E  +  00 

2.27  90E  +  03 

2 

.  0000E  +  02 

112 

2 . 3000E  +  20  -2 . 17  00E  +  00 

2 . 1330E  +  03 

2 

.  5000E  +  02 

113 

9 . 1100E  +  1 9  -2 . 0400E  +  00 

2 . 01 60E+03 

3 

.  0000E  +  02 

114 

4 . 2200E+19  -1 . 9400E+00 

1 . 9170E+03 

3 

.  5000E  +  02 

115 

2 . 1 900E+1 9  -1 . 8500E+00 

1 . 832 0E  +  03 

4 

.  0000E  +  02 

116 

1 . 2400E  +  19  -1 . 77  00E  +  00 

1 . 7  5  90E  +  03 

4 

.  5000E  +  02 

117 

7 . 5500E  +  18  -1 . 7  000E  +  00 

1 . 6940E  +  03 

5 

.  0000E  +  02 

118 

3 . 2500E  +  18  -1 . 5800E  +  00 

1 . 5830E  +  03 

6 

.  0000E  +  02 

119 

1 . 6200E  +  18  -1 . 4 900E  +  00 

1 . 4 92 0E  +  03 

7 

.  0000E  +  02 

120 

9 . 0500E  +  17  -1 . 4100E  +  00 

1 . 4140E  +  03 

8 

.  0000E  +  02 

121 

5 . 4700E  +  17  -1 . 3400E  +  00 

1.347  0E  +  03 

9 

.  0000E  +  02 

122 

3 . 5300E  +  17  -1 . 2800E  +  00 

1.2880E+03 

1 

.  0000E  +  03 

123 

7 . 1400E  +  16  -1 . 0 600E  +  00 

1 . 07  30E  +  03 

1 

.  5000E  +  03 

124 

2 . 5300E  +  1 6  -9 . 1000E-01 

9 . 3200E  +  02 

2 

.  0000E  +  03 

43. 

N2H2+NH2=NH3+NNH 

1. 

000E  +  13 

0.00 

1000.0 

ee 

44  . 

NH3+0H=NH2+H20 

2  . 

040E  +  06 

2.04 

566.0 

ee 

45. 

NH3 (+M) =NH2+H (+M) 

6. 

600E  +  17 

0.00 

97280.0 

ee 

Low  pressure  limit:  0.20400E+16 

0 . 00000E  +  00 

0 . 78780E+05 

46. 

HN0+0H=N0+H20 

1. 

295E+07 

1.88 

-958.0 

ee 

47  . 

HNO+NH2=NH3+NO 

2  . 

000E  +  13 

0.00 

1000.0 

ee 

48. 

NNH (+M) =N2+H (+M) 

4  . 

100E  +  09 

1.13 

5186.0 

ee 

Low  pressure  limit:  0.10000E+14 

0 . 50000E+00 

0 . 30600E+04 

N20  Enhanced  by 

5 . 000E  +  00 

H20  Enhanced  by 

9 . 000E  +  00 

N2 

Enhanced  by 

1 . 000E  +  00 

02 

Enhanced  by 

8 . 200E-01 

NH3  Enhanced  by 

5 . 000E  +  00 

49. 

HN0+02=N0+H02 

1 . 

373E  +  01 

3.56 

7932.3 

ff 

50. 

HOOCL= 

=HCL+02 

8. 

602E  +  12 

0.15 

34.0 

gg 

51. 

HOOCL= 

=HOCLO 

2  . 

659E  +  12 

0.35 

39.8 

gg 

52  . 

0H+CL02 (+M) =H0CL02 (+M) 

1. 

000E  +  00 

0.00 

0.0 

hh 

105 

1 . 2 600E  +  23  -3 . 9400E  +  00 

5 . 4100E  +  02 

1 

.  0000E  +  01 

106 

2.5800E+23  -3.9200E+00 

6 . 7  600E  +  02 

2 

.  5000E  +  01 

107 

4 . 1200E  +  23  -3 . 8 900E  +  00 

8 . 0500E  +  02 

5 

.  0000E  +  01 

108 

5 . 2500E  +  23  -3.8800E  +  00 

8 . 9100E+02 

7 

.  5000E  +  01 

109 

6 . 1500E  +  23  -3 . 8 600E  +  00 

9 . 57  00E  +  02 

1 

.  0000E  +  02 

110 

7.5000E+23  -3 . 8400E+00 

1 . 05  60E  +  03 

1 

.  5000E  +  02 

Ill 

8 . 4 600E  +  23  -3.8200E  +  00 

1 . 1290E  +  03 

2 

.  0000E  +  02 

112 

9 . 1 600E  +  23  -3 . 8000E  +  00 

1 . 1880E+03 

2 

.  5000E  +  02 

113 

9 . 6800E+23  -3.7800E+00 

1 . 237  0E  +  03 

3 

.  0000E  +  02 

114 

1 . 0100E  +  24  -3 . 77  00E  +  00 

1 . 27  80E  +  03 

3 

.  5000E  +  02 

115 

1 . 0300E  +  24  -3 . 7 600E  +  00 

1 . 3150E+03 

4 

.  0000E  +  02 

116 

1 . 0500E  +  24  -3 . 7500E  +  00 

1.347  0E  +  03 

4 

.  5000E  +  02 

117 

1 . 0700E  +  24  -3 . 7300E  +  00 

1.37  60E  +  03 

5 

.  0000E  +  02 

118 

1 . 0800E  +  24  -3 . 7100E  +  00 

1 . 4250E  +  03 

6 

.  0000E  +  02 

119 

1 . 0800E  +  24  -3 . 6900E  +  00 

1 . 4 67 0E  +  03 

7 

.  0000E  +  02 

120 

1 . 0700E  +  24  -3 . 6800E  +  00 

1 . 5030E  +  03 

8 

.  0000E  +  02 

Approved  for  public  release;  distribution  is  unlimited. 


34 


121 

1 . 0500E+24  - 

3. 6600E+00 

1 . 5350E+03 

9 . 0000E+02 

122 

1.0300E+24  - 

3. 6500E+00 

1 . 5  630E  +  03 

1.0000E+03 

123 

8 . 9900E+23  - 

3 . 5800E+00 

1 . 6660E+03 

1 . 5000E+03 

124 

7 . 6700E+23  - 

3 . 5200E+00 

1 . 7330E+03 

2 . 0000E  +  03 

53. 

CL0+CL0=CL2+02 

6 . 560E  +  10 

0.66 

3759.4 

ii 

54  . 

CLO+CLO=CL+CLOO 

8.190E+10 

0.77 

4307.8 

ii 

55. 

CL0+CL0=CL02+CL 

3.77  OE  +  13 

0.01 

5754 . 4 

ii 

56. 

CLOO (+M) =CL+02 (+M) 

1 . 000E+00 

0.00 

0.0 

j  j 

105 

6 . 2 600E  +  0 9 

8 . 0000E-02 

2 . 5080E+03 

1 . 0000E  +  01 

106 

1 . 5600E  +  10 

8 . 0000E-02 

2 . 512  0E  +  03 

2 . 5000E  +  01 

107 

3 . 1300E  +  10 

8 . 0000E-02 

2 . 51 90E+03 

5 . 0000E  +  01 

108 

4 . 6700E+10 

8 . 0000E-02 

2 . 52  60E  +  03 

7 . 5000E  +  01 

109 

6 . 2100E  +  10 

9 . 0000E-02 

2 . 532  0E  +  03 

1 . 0000E  +  02 

110 

9 . 2300E  +  10 

9 . 0000E-02 

2 . 5430E+03 

1 . 5000E  +  02 

Ill 

1 . 2200E  +  11 

9 . 0000E-02 

2 . 5530E+03 

2 . 0000E  +  02 

112 

1 . 5100E  +  11 

9 . 0000E-02 

2 . 5  630E  +  03 

2 . 5000E  +  02 

113 

1 . 7  800E  +  11 

9 . 0000E-02 

2.572  0E  +  03 

3 . 0000E  +  02 

114 

2 . 0600E  +  11 

9 . 0000E-02 

2 . 5810E+03 

3 . 5000E  +  02 

115 

2 . 3200E  +  11 

9 . 0000E-02 

2 . 5880E+03 

4 . 0000E  +  02 

116 

2 . 6800E+11 

9 . 0000E-02 

2 . 6040E+03 

4 . 5000E  +  02 

117 

2 . 8200E  +  11 

1 . 0000E-01 

2 . 6020E+03 

5 . 0000E  +  02 

118 

3 . 3200E  +  11 

1 . 0000E-01 

2 . 6160E+03 

6 . 0000E  +  02 

119 

3 . 7500E  +  11 

1.0000E-01 

2 . 6270E+03 

7 . 0000E  +  02 

120 

3 . 8500E  +  11 

1 . 2000E-01 

2 . 6260E+03 

8 . 0000E  +  02 

121 

4 . 6200E  +  11 

1.1000E-01 

2 . 6490E+03 

9 . 0000E  +  02 

122 

4 . 9700E+11 

1 . 1000E-01 

2 . 6570E+03 

1 . 0000E  +  03 

123 

6 . 8100E  +  11 

1 . 2000E-01 

2 . 6990E+03 

1 . 5000E  +  03 

124 

8 . 2  900E  +  11 

1 . 3000E-01 

2 . 7300E+03 

2 . 0000E  +  03 

57  . 

CL+NH3=NH2+HCL 

3.24  8E  +  05 

2.53 

1322.0 

kk 

58. 

CLO+NH3=NH2+CLOH 

1 . 130E  +  00 

3.85 

8631.5 

11 

59. 

CL03+NH3=NH2+H0CL02 

8.190E+09 

1.01 

4480.7 

11 

60. 

CLO+NH2=HCL+HNO 

2 . 830E  +  16 

-1.08 

256.3 

mm 

61. 

CL0+NH2=CL+NH20 

1.020E+15 

-0.62 

47 . 7 

mm 

62. 

CL03+NH2=CL02+NH20 

5. 960E  +  15 

-0 .47 

47 . 7 

mm 

63. 

CLO+NH20=CLOH+HNO 

2 . 000E  +  13 

0.00 

0.0 

nn 

64  . 

CL02+NH20=H0CL0+HN0 

2 . 000E  +  13 

0.00 

4210 . 0 

oo 

65. 

CL03+NH20=H0CL02+HN0 

2 . 000E  +  13 

0.00 

0.0 

nn 

6  6 . 

OH+CL03 (+M) =H02+CL02 (+M) 

1 . 000E  +  00 

0.00 

0.0 

PP 

105 

1 . 2  900E  +  14 

9 . 0000E-02 

3. 9000E+01 

1.0000E+01 

124 

1 . 2 900E+14 

9 . 0000E-02 

3. 9000E+01 

2 . 0000E  +  03 

67  . 

CLO+NO (+M) =CL+N02 (+M) 

1 . 000E  +  00 

0.00 

0.0 

qq 

105 

1 . 3400E  +  13  - 

5 . 0000E-02 

1.8000E+01 

1.0000E+01 

124 

1 . 3400E  +  13  - 

5 . 0000E-02 

1 . 8000E+01 

2 . 0000E  +  03 

68. 

CL02+N0=CL0+N02 

3 . 120E  +  11 

0.39 

-761.0 

rr 

69. 

NOCL+H20=HONO+HCL 

1 . 224E  +  11 

0.23 

11610.0 

ss 

70. 

NOCL+NO=CLOO+N2 

1 . 7  00E  +  13 

0.00 

5625.0 

bb 

71. 

HCL+CL=CL2+H 

1 . 500E+15 

0.00 

47550.0 

tt 

72  . 

CLOH+HNO=HCL+HONO 

3 . 480E  +  02 

3.06 

6068.3 

uu 

73. 

HCL04= 

=CL03+OH 

1 . 500E+17 

0.00 

52655.5 

PP 

74  . 

NH4CL04=0 JHNH3+CL03 

1 . 500E  +  17 

0.00 

61560.0 

vv 

75. 

0JHNH3=NH2+H20 

1.091E+16 

0.47 

7036.9 

ss 

76. 

OH+NH3=0 JHNH3 

7 . 286E  +  10 

0.64 

-865.1 

bb 

77  . 

HCL04+H0N0=H20+0N0CL03 

3.226E+10 

0.19 

-1.0 

ss 

78. 

0N0CL03=CL02+N0+02 

2 . 458E  +  12 

0.35 

17353.0 

ss 

79. 

NH2+NO=N2+OH+H 

8 . 990E+11 

0.00 

0.0 

ee 

80. 

N2H3+N2H3=N2H4+N2H2 

1 . 200E+13 

0.00 

0.0 

ww 

81. 

N2H4+N02=N2H3+H0N0 

5. 603E+04 

2.54 

7787 . 0 

ss 

82  . 

N2H4+CL03=N2H3+H0CL02 

8.190E+09 

1.01 

0.0 

bb 

83. 

NH20+OH (+M) =NH2+H02 (+M) 

1 . 000E  +  00 

0.00 

0.0 

XX 

105 

2 . 0000E  +  24  - 

2 . 9900E+00 

2 . 82  96E  +  04 

1.0000E+01 

123 

2 . 0000E  +  24  - 

2 . 9900E+00 

2 . 8296E+04 

1 . 5000E  +  03 

124 

1 . 9900E+24  - 

2 . 9900E+00 

2 . 82  96E  +  04 

2 . 0000E  +  03 

84  . 

NH2+N02=N0+NH20 

9 . 030E  +  11 

0.03 

-1512 . 1 

YY 
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85.  NH2+N02 (+M) =H2NN02 (+M)  l.OOOE+OO  0.00  0.0 

105:  2 . 8 900E  +  2 9  -5.6300E  +  00  3.5720E  +  03  1.0000E  +  01 

106:  1 . 5000E  +  2 9  -5.4400E+00  3.5050E  +  03  2.5000E  +  01 

107:  5 . 3800E+28  -5.2300E+00  3.4070E+03  5.0000E+01 

108:  2 . 4700E+28  -5.0800E+00  3.3390E+03  7.5000E+01 

109:  1 . 3100E+28  -4.9700E+00  3.2850E+03  1.0000E+02 

110:  4 . 7200E+27  -4.7900E+00  3.1990E+03  1.5000E+02 

Ill:  2 . 0700E+27  -4.6600E+00  3.1300E+03  2.0000E+02 

112:  1 . 0300E+27  -4.5500E+00  3.0690E+03  2.5000E+02 

113:  5 . 6200E+2 6  -4.4500E+00  3.0160E  +  03  3.0000E  +  02 

114:  3 . 27  00E  +  2 6  -4.3600E+00  2.9670E+03  3.5000E  +  02 

115:  2 . 0000E  +  2 6  -4.2900E+00  2.9230E  +  03  4.0000E  +  02 

116:  1 . 2800E  +  2 6  -4.2200E  +  00  2.8810E  +  03  4.5000E  +  02 

117:  8 . 5000E+25  -4.1600E+00  2.8430E+03  5.0000E+02 

118:  4 . 0700E+25  -4.0400E+00  2.7730E+03  6.0000E+02 

119:  2 . 1300E+25  -3.9500E+00  2.7100E+03  7.0000E+02 

120:  1 . 2000E+25  -3.8600E+00  2.6540E+03  8.0000E+02 

121:  7 . 1300E+24  -3.7800E+00  2.6020E+03  9.0000E+02 

122:  4 . 4400E+24  -3.7100E+00  2.5540E+03  1.0000E+03 

123:  6 . 6700E+23  -3.4300E+00  2.3580E+03  1.5000E+03 

124:  1.6500E+23  -3.2300E+00  2.2090E+03  2.0000E+03 


86. 

H2NN02+NH2=HN JN02+NH3 

4 . 820E  +  00 

3.60 

770 . 0 

aaa 

87  . 

HN JN02+CL02=HN0 JN02+CLO 

1 . 000E  +  13 

0.00 

0.0 

bb 

00 

00 

HNO JN02=HN0+N02 

3.75  6E  +  13 

0.31 

12343.0 

ss 

89. 

NH20+N02=HN0+H0N0 

5. 603E+04 

2.54 

0.0 

ss 

90. 

CL02+HNO=HOCLO+NO 

1 . 2  95E  +  01 

3.24 

651.5 

ss 

91. 

H0CL02+CL=CL03+HCL 

3 . 550E-01 

4.07 

-340.0 

bbb 

92. 

0+0H+M=H02+M 

5 . 000E  +  1 6 

0.00 

0.0 

ccc 

93. 

NH3+NH3=DINH3 

3.292E+07 

1.82 

-1391.5 

ddd 

94  . 

NH4CL04+NH3=AP  NH3 

3.292E+07 

1.82 

-1391.5 

ddd 

95. 

DINH3+HCL04=AP_NH3 

3.292E+07 

1.82 

-1391.5 

ddd 

96. 

C0+02=C02+0 

2 . 530E+12 

0.00 

47688.0 

ee 

97  . 

0H+CH20<=>CH0+H20 

3.430E+09 

1.20 

-447.0 

ee 

98. 

C2H3+02<=>CH0+CH20 

4.580E+16 

0 

5 - 1 

1 

1015.0 

ee 

99. 

CHO+M=H+CO+M 

1 . 870E+17 

O 

O 

5 - 1 

1 

17000.0 

fff 

100. 

C^CCL^O+OH=VCDJO+CLOH 

3 . 500E  +  12 

0.00 

22810.0 

eee 

NOTE:  A  units  mole-cm-sec-K,  E  units  cal/mole 


Notes 


a.  Estimate  obtained  via  the  Quantum-Rice-Ramsperger-Kassel  (QRRK)  method:  5-ethenyl- 
2,8-decadiene  — >  l-butene-3yl  +  2,6-octadiene-lyl  (Chen  and  McQuaid  2009). 

b.  Estimate  obtained  via  QRRK:  CHsCHCHOEOfcCHCHC'Hi  CH3CHCHC*H2  + 
CH2CHCHCH2. 

c.  Estimated  from  trans-C H 2C' H C H C' H C' HC H 2  — >  1,3-cyclohexadiene  (Grimme  et  al.  1981). 

d.  From  (Dean  and  Bozzelli  2000). 

e.  Estimate  based  on  rate  expression  for  C2H2+C2H3  — >  CtCCHCHC*!!  found  in  Benson 
(1989). 

f.  From  Kiefer  et  al.  (1988). 

g.  Estimate  based  on  transition  state  theory  (TST).  BIK-determined  entropies  were  employed  to 
calculate  A.  G3//BlK-based  results  were  employed  to  calculate  Ea. 

h.  From  Weissman  and  Benson  (1988). 

i.  From  Tsang  and  Herron  (1991). 

j.  Estimated  from  a  rate  expression  for  C3H6  +  Cl  — »  CFFCHC^Fh  +  HC1  obtained  from  Pilgrim 
and  Taatjes  (1997). 
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k.  Estimated  from  a  TST-derived  rate  expression  for  CIO  +  c/ .s’  -  C  H  3  C  H  ( C  2  H  3 )  C  H  3  — »  ClOH  + 
CH3C*(C2H3)CH3.  Density  functional  theory  (DFT)-determined  entropies  were  employed  to 
calculate  A.  CBS-QB3-based  results  were  employed  to  calculate  Ea. 

l.  Estimated  from  a  rate  expression  for  C2H4  +  Cl  — >  C*H2CH2C1  obtained  from  Knyazev  et  al. 
(1999). 

m.  Estimated  from  a  rate  expression  for  C3H6  +  Cl  — >  CFEOHCFhCl  obtained  from  Kaiser  and 
Wallington  (1996). 

n.  Estimated  from  a  TST-derived  rate  expression  for  CH3CH(CHCH2)CH2C*HCHC1CH3  -> 
CH3CHCHC»H2  +  CH2CHCHC1CH3.  DFT-determined  entropies  were  employed  to  calculate 
A.  CBS-QB3-based  results  were  employed  to  calculate  Ea. 

0.  Estimated  from  a  TST-derived  rate  expression  for  CH3CHCHCH2CH(CH3)C»HCH2C1  — > 
CH3CHCHC»H2  +  CH3CHCHCH2C1.  DFT-determined  entropies  were  employed  to  calculate 
A.  CBS-QB3-based  results  were  employed  to  calculate  Ea. 

p.  Estimated  via  QRRK  with  a  variational  transition  state  theory  (VTST)-derived  rate 
expression  for  CH2CHC*HCH3  +  C102  ->  CH2CHC0C10HCH3  (5.94  x  104 

T2. 367 ' 7exp(!  723.4/^7)  cm3  mole"1  s'1)  and  a  TST-derived  rate  expression  for 
CH2CHC0C10HCH3  ->  CH2CHCOCH3  +  ClOH  (6.63  x  1011  7<I21524exp(-8639.3/R7)  s"1). 

q.  From  Curran  et  al.  (2004). 

r.  Estimated  based  on  a  TST-derived  rate  expression  for  CH3CHCHCH(CHCH2)CH2C*H2  — » 
C2HCH2+  CH3CHCHC*HCHCH2.  BIK-determined  entropies  were  employed  to  calculate  A. 
G3MP2//BlK-based  results  were  employed  to  calculate  Ea. 

s.  Estimated  based  on  microscopic  reversibility  (MR)  and  a  reverse  rate  expression  for  CH3  + 
C2H4  — >  CH3CHC*H2  obtained  from  Curran  (2006). 

t.  Estimated  based  on  a  QRRK-derived  rate  expression  for  CH2CHC»H2  +  N02  — » 
CH2CHCH20*  +  NO.  Input  included  a  VTST-derived  rate  expression  for  CH2CHC*H2  +  N02 
— »  CH2CHCH2ONO  (1.87  x  109  7°  95983exp(-948/i?7)  cm3mole'1s'1)  and  a  rate  expression  for 
CH2CHCH2ONO  — >  CH2CHCH20»  +  NO  that  was  estimated  based  on  MR  and  a  reverse  rate 
expression  for  CH30  +  NO  — »  CH3ONO  obtained  from  Caralp  et  al.  (1998). 

u.  Estimated  based  on  a  TST-derived  rate  expression  for  CH2CHCO*HCH2CHCH2  — » 
CH2CHCHO  +  CH2CHC*H2.  BIK-determined  entropies  were  employed  to  calculate  A. 
G3//BlK-based  results  were  employed  to  calculate  Ea. 

v.  Estimated  based  on  a  TST-derived  rate  expression  for  CH2CHCO*HCH2CH2CH3  — » 
CH2CHCHO  +  CH3CH2C»H2.  BIK-determined  entropies  were  employed  to  calculated. 
G3//BlK-based  results  were  employed  to  calculate  Ea. 

w.  Estimated  based  on  a  QRRK-derived  rate  expression  for  CH3C*H2  +  N02  — »  CH3CH20»  + 
NO.  Input  included  a  rate  expression  for  CH3C*H2  +  N02  — >  CH3CH2ONO  that  was 
estimated  to  be  2.03  x  1013  cm3mol'1s'1  and  a  rate  expression  for  CH3CH2ONO  — »  CH3CH20* 
+  NO  that  was  estimated  based  on  MR  and  a  reverse  rate  expression  for  CH3O  +  NO  — > 
CH3ONO  obtained  from  Caralp  et  al.  (1998). 

x.  Estimated  based  on  a  TST-derived  rate  expression  for  CH2CHCH2CH20*  — »  CH2CHC*H2  + 
CH20.  BIK-determined  entropies  were  employed  to  calculate  A.  G3//BlK-based  results 
were  employed  to  calculate  Ea. 

y.  Estimated  based  on  a  TST-derived  rate  expression  for  CH3CHCHCH2CH2CHCHC*H2  — > 
CH3CHCHC»H2  +  CH2CHCHCH2.  BIK-determined  entropies  were  employed  to  calculate  A. 
G3//BlK-based  results  were  employed  to  calculate  Ea. 
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z.  From  Tsang  and  Hampson  (1986). 

aa.  Estimated  from  a  rate  expression  for  C2H5  +  O  — »  C2H4  +  OH  obtained  from  Harding  et  al. 
(2005). 

bb.  Estimated  in  this  study. 

cc.  Estimated  from  a  TST-derived  rate  expression  for  CH2CHCH2CHO-CH3  CH2CHCHO 
+  CH3.  BIK-determined  entropies  were  employed  to  calculate  A.  G3//BlK-based  results 
were  employed  to  calculate  Ea. 

dd.  Estimated  based  on  a  TST-derived  rate  expression  for  CH2CHCH2CHO»CH3  — > 

CH2CHCOCH3  +  H.  BIK-determined  entropies  were  employed  to  calculate  A.  G3//B1K- 
based  results  were  employed  to  calculate  Ea. 
ee.  From  Anderson  et  al.  (2011). 

ff.  Estimated  via  TST  with  an  H-atom  tunneling  correction.  DFT-determined  entropies  were 
employed  to  calculate  A.  G4-  and  CBS-QB3-based  results  were  employed  to  calculate  Ea. 
gg.  From  Xu  et  al.  (2003) 

hh.  Estimated  via  QRRK  with  input  obtained  from  Xu  et  al.  (2003). 
ii.  From  Zhu  and  Lin  (2003a). 

jj.  Estimated  via  QRRK  with  input  obtained  from  Zhu  and  Lin  (2003b). 
kk.  From  Monge-Palacios  and  Espinosa-Garcia  (20 1 0). 

II.  From  Xu  and  Lin  (2007) 
mm.  From  Zhu  and  Lin  (2007). 

nn.  Estimated  from  results  obtained  for  Cl  +  NH20  — »  HC1  +  HNO. 

00.  A  was  estimated  based  on  results  obtained  for  Cl  +  NH20  — >  HC1  +  HNO.  G4-based  results 
were  employed  to  calculate  Ea. 

pp.  Estimated  via  QRRK  with  input  obtained  from  Zhu  and  Lin  (2001). 
qq.  From  Zhu  and  Lin  (2004). 

rr.  Estimated  based  on  a  rate  expression  for  CIO  +  NO  — >  Cl  +  N02  obtained  from  Zhu  and  Lin 
(2004). 

ss.  Estimated  based  on  TST  with  an  H-atom  tunneling  correction.  DFT-determined  entropies 
were  employed  to  calculate  A.  G4-based  results  were  employed  to  calculate  Ea. 
tt.  NIST  (Manion  et  al.  2013)  Two-parameter  fit  to  rates  for  HC1  +  Cl  — >  C12  +  H. 
uu.  From  Xu  and  Lin  (2010). 

vv.  A  was  estimated  based  on  the  rate  expression  for  HC104  — >  C103  +  OH.  G4-based  results 
were  employed  to  calculate  Ea. 
ww.  From  Sun  et  al.  (2009). 

xx.  Estimated  via  QRRK  with  a  rate  expression  for  NH20  +  OH  — >  NH200H  that  was 

estimated  to  be  1.81  x  1013  cm3mole-ls-l.  Input  included  a  rate  expression  for  NH200H 
— >  NH2  +  H02  that  was  based  on  MR  and  reverse  rate  expression  for  NH2  +  H02  — > 
NH200H  that  was  obtained  from  Bozzelli  and  Dean  (1989). 
yy.  From  Klippenstein  et  al.  (2013). 

zz.  Estimated  via  QRRK  with  input  obtained  from  Klippenstein  et  al.  (2013). 

aaa.  Estimated  based  on  a  rate  expression  for  N2H4  +  NH2  — >  NH3  +  N2H3  obtained  from  Li  and 
Zhang  (2006). 

bbb.  Estimated  based  on  a  rate  expression  for  HOC1  +  Cl  — >  HC1  +  CIO  obtained  from  Wang  et 
al.  (2003). 

ccc.  From  Burke  et  al.  (2010). 
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ddd.  Estimated  based  on  a  rate  expression  for  NH3  +  HC104  — »  NH4C104  obtained  from  Zhu 
and  Lin  (2008). 

eee.  From  Ho  et  al.  (1995)  and  Procaccini  et  al.  (2000). 

fff.  From  Timonen  et  al.  (1987). 
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Appendix  B.  Coefficients  for  Calculating  Thermodynamic 
Property  Estimates  for  Species  in  Mechanism  R2 
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Coefficients  for  computing  estimates  for  the  heat  capacity  [cp(T)],  enthalpy  H(T),  and  entropy 
[S(T)]  of  species  in  mechanism  R2  are  given  in  Table  B-l.  The  data  conform  to  standard  Chemkin 
II  formatting  rules.1  Because  the  CHEMKIN  II  framework  limits  the  labels  for  species  to  16  or 
fewer  characters,  for  many  species  it  was  not  possible  to  assign  labels  that  can  be  unambiguously 
translated  into  a  molecular  structure.  A  look-up  table  that  performs  that  function  is  provided  in 
Appendix  C. 


Table  B-l  Coefficients  for  calculating  thermodynamic  property  estimates  for  species  in  mechanism  R2 


THERMO 

300.000  1500.000  5000.000 

CHO  ESTC  1H 

3. 69472521E+00  3 . 185942 96E-03-1 
3 . 82240388E  +  03  4 . 69145660E+00  3 
-1 . 72919680E-09  3 . 98120351E-13  4 
CH20  L  8/88H  2C 

1 . 7 60 69008E  +  00  9 . 20000082E-03-4 
-1.39958323E+04  1 . 36563230E  +  01  4 
-3.79285261E-08  1 . 31772 652E-11-1 
CH3  121286C  1H 

0 . 02844052E  +  02  0 . 06137974E-01-0 
0 . 01643781E  +  06  0 . 05452 697E+02  0 
0 . 01621829E-06-0 . 05864953E-10  0 
CO  121286C  10 

0 . 03025078E  +  02  0 . 01442689E-01-0 
-0 . 01426835E  +  06  0 . 06108218E  +  02  0 
0 . 05581 94 4E-07-0 . 02474 951E-1 0-0 
C02  121286C  10 

0 . 04453623E  +  02  0 . 03140169E-01-0 
-0. 0489669 6E+06-0. 0955395 9E+01  0 
0 . 06866687E-07-0 . 02117280E-10-0 
C2H2  121386C  2H 

0 . 04436770E  +  02  0 . 0537 6039E-01-0 
0 . 025667 6 6E +  06-0 . 02800338E  +  02  0 
0 . 09078992E-07-0 . 01912746E-10  0 
C2H3  12787C  2H 

0 . 05933468E  +  02  0 . 0401774 6E-01-0 
0 . 03185435E  +  0  6-0 . 08530313E  +  02  0 
-0 . 01321642E-07-0 . 01184784E-10  0 
C*CCL*0  6/19/97  C  3H 

1 . 14618122E  +  01  8. 4510946 9E -03-2 
-2 . 42075137E  +  04-3 . 15944471E  +  01  2 
1 .23566507E-08-2 . 2 072384 OE-12-2 
VCD JO  VCD JO_G3MP2_G3C  3H 

8 . 7 6612341E  +  00  7 . 77928992E-03-2 
7 . 7990855 6E +03-1 . 86491247E+01  2 
1.29707093E-08-2.52092364E-12  9 
VCDO  VCDO_G3MP2_G3_C  3H 

9 . 39918255E  +  00  1 . 01168249E-02-3 
-1 .23594302E  +  04-2 . 47300373E  +  01  1 
6. 98466320E-09-9. 7 9237748E-13-9 


10  1  0G 
08841412E-06 
53025733E+00 
08521632E+03 
10  1  00G 
42258813E-0  6 
7  9372315E  +  00- 
43089567E+04 
3  G 

02230345E-04 
02430443E+02 
01 64237  8E  +  0  6 

1  G 
05630828E-05 
032  62452E  +  02 
01431054E  +  0  6 

2  G 
01278411E-04 
02275725E+02 
04837314E+06 

2  G 
01912817E-04 
02013562E+02 
02  612  4  4  4E  +  0  6 

3  G 

03966740E-05- 
0245927  6E  +  02 
03335225E  +  0  6 
30  1CL  1G 

9007387  9E-0  6 
31820028E+00 
1301 6608E+04 
3N  OO  1G 

4  902  67  90E-0  6 

10804239E+00 
65945852E+03 
4N  OO  1G 

36959545E-06 
13047  84  6E  +  00 
4  94  90843E  +  03 


300.000  5000 

1 . 68761454E-10- 
1.88364239E-03 
6 . 234  92345E  +  00 
300.000  3500 

1.00641212E-09- 
9. 90833369E-03 
6 . 02812  900E-01 
0300.00  5000 

0.03785161E-08- 
0 . 01112410E+00- 
0 . 06789794E+02 
0300.00  5000 

0.01018581E-08- 
0 . 01511941E-01- 
0 . 04848897E+02 
0300.00  5000 

0.02393997E-08- 
0 . 09922072E-01- 
0 . 0101884  9E+03 
0300.00  5000 

0.03286379E-08- 
0 . 01519045E+00- 
0 . 08805378E+02 
0300.00  5000 

0 . 014412  67E-08 
0.0737147  6E-01 
0 . 01155620E+03 
300.000  5000 

4.51447285E-10- 
3.25587682E-02- 
1 . 654  6502  9E  +  01 
300.000  5000 

3. 69252827E-10- 
2.79698291E-02- 
1 . 54835048E+01 
300.000  5000 

5. 1417294 6E -10- 
2 . 922  6987  9E-02- 
1 . 97052158E+01 


.000  1367.000 
9.77  966305E-15 
1 . 78452098E-06 

.000  1000.000 
8 . 83855640E-14 
3 . 73220008E-05 
1 . 00197170E+04 
.00  1000.00 
0 . 02452159E-12 
0 . 01680220E-03 

.00  1000.00 
0 . 06910952E-13 
0 . 03881755E-04 

.00  1000.00 
0 . 01 669033E-12 
0 . 01040911E-03 

.00  1000.00 
0 . 02156710E-12 
0 . 01 61 631 9E-03 

.00  1000.00 
0 . 02378644E-12 
0 . 02109873E-04 

.000  1398.000 
2 . 6237  6318E-14 
2 . 78476237E-05 

.000  1403.000 
2 . 07329959E-14 
2 . 62572868E-05 

.000  1392.000 
2 . 94835284E-14 
1 . 99974542E-05 


01 

2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 

11 

2 

3 

4 

01 

2 

3 

4 

01 

2 

3 

4 


1  Kee  RJ,  Rupley  FM,  Miller  JA.  Chemkin  II:  A  Fortran  chemical  kinetics  package  for  the  analysis  of  gas-phase  chemical 
kinetics.  Albuquerque  (NM):  Sandia  National  Laboratories;  1989.  Report  No.:  SAND89-8009. 
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CTCV  10/  2/9  C  4H  4  0  0G 

1.123677  61E+01  8 . 823  64166E-03-3 . 071532  01E-0  6 
2 . 92566868E+04-3. 5733636 6E +01-2 . 9996813 6E-02 
1 . 423224 91E-08-2 . 5032395 6E -12  3 . 28915781E+04 
CDCCDCJ  CDCCDC_BK1+VINC  4H  5N  OO  OG 

1.131 635 93E+01  1 . 05404088E-02-3 . 49684912E-06 
3 . 77210212E+04-3 . 5 9442 045E+01-2 . 23870 963E+00 
2 . 132 6873 6E-08-3. 90283982E-12  4 . 16788073E+04 
CDCCDC  CDCCDC_BK1 . outC  4H  6N  OO  OG 

1 . 14271967E+01  1 . 2 95844 95E-02-4 . 32307232E-06 
7 . 61440843E+03-3 . 84271070E+01-3 . 16530 04 8E+00 
2 . 04344500E-08-3. 60678263E-12  1 . 20575132E+04 
CCDOV  9/20/  9  THERMO  4H  60  1  OG 

1 . 11282775E+01  1 . 59830571E-02-5 . 50607262E-06 
-1 . 91689528E+04-3.24493094E+01  2 . 5457 6818E+00 
-1 . 18 684 688E-0 9  9 . 30  612  7  08E-13-1 . 550  6814  9E  +  04 
CDCCJC  CDCC JC_BK1 . outC  4H  7N  OO  OG 

1 . 01 66114 1E+01  1 . 66497350E-02-5 . 60871682E-06 
1 . 06328953E+04-2 . 9533747 6E+01  1 . 5030 634 4E-02 
8.25242412E-09-1 . 13970812E-12  1 . 4284254 9E+0 4 
VCOJC  VCO JC_G3MP2_BKC  4H  7N  OO  1G 

1 .26649548E+01  1 . 66701855E-02-5 . 62013651E-06 
1.02185455E+03-3. 98930857E+01  1 . 0018434 6E+00 
1 .24303453E-08-2 . 01042 64 4E-12  5 . 00471004E+03 
CJW  CW+C*CC*CCJ  C  5H  7N  OO  0G 

1 . 42616569E+01  1 . 598014 97E-02-5 . 53201590E-06 
1. 60483386E+04-5.14843851E+01-3.16510508E+00 
1 . 7318114 6E-08-2 . 70632680E-12  2 . 1 90 65453E+04 
WCOJ  10/  1/  9  THERMO  5H  70  1  OG 

1 . 6158  98  92E  +  01  1 . 62  62  8  904E-02-5 . 4 9172  54  3E-0 6 
1 . 04839262E+04-5. 7728066 9E +01-4 . 66489790E-01 
2 . 04451484E-08-3 . 53334 092E-12  1 . 58832224E+04 
CW  CW_BK1  .  out  C  5H  8N  OO  0G 

1 . 38617204E+01  1 . 7 9152384E-02-5 . 990177 88E-0 6 
2 .33368455E+03-4 . 95140504E+01-2 . 58120389E+00 
1 . 85895355E-08-3 . 00239718E-12  7 . 64306740E+03 
CJCVC  E_C JCVC_BK1 . ouC  5H  9N  OO  0G 

1 .263117 14E+01  2 . 11 968230E-02-7 . 14277909E-06 
1 . 47076618E+04-3.78871255E+01  2 . 63277 613E+00 
4 . 0230732 9E-0 9-1 . 64144521E-13  1 . 8 6551000E+04 
CVCCVCL  GA  C  8H  13  0CL  1G 

2 . 4  6598556E  +  01  3 . 103  687  4  4E-02-1 . 0  65  667  8  6E-05 
-1 .22695766E+04-9. 68 68830 6E+0 1-2 . 27 668100E+00 
2 . 83739559E-08-4 . 6302 182 8E-12-3 . 08132 966E+03 
VY6DE13J  10/  1/  9  THERMO  8H  9  0  0G 

2 . 16368750E+01  2 . 26941806E-02-7 . 71538060E-06 
1 . 9655938 6E +04-9. 23086148E+01-7 . 4 93602 60E+00 
4 . 54442 654E-08-8 . 4 1355134E-12  2 . 83839413E+04 
WVJV  10/  1/9  C  8H  9  0  0G 

2 . 2  6301412E  +  01  2 . 15222 007E-02-7 . 42  0  962  8  6E-0  6 
4 . 14479674E+04-9. 30067034E+01-2 . 38964234E+00 
4 . 44147356E-08-8 . 62995317E-12  4 . 88 90438 9E+04 
VWIV  9/24/  9  THERMO  8H  10  0  0G 

2 . 3522312  9E  +  01  2 . 263  9130  6E-02-7 . 7507  84  45E-0  6 
1 . 48137893E+04-9. 98670223E+01-4 . 64836883E+00 
5.31747258E-08-1 . 04157863E-11  2 . 28979384E+04 
VY6DE13  10/  1/  9  THERMO  8H  10  0  0G 

2 . 1717  97  98E  +  01  2 . 52  4  72190E-02-8 . 7  04  62  98  9E-0  6 
8 . 6754082 6E+03- 9 . 4 9443753E+01-7 . 2 64 617 99E+00 
3. 96451 662E-08-7 . 10924565E-12  1 . 78323550E+04 


300.000  5000.000  1392.000 

4 . 82681345E-10-2 . 82467653E-14 
3 . 79242319E-02-3 . 25061220E-05 
2 . 37  91222  9E  +  01 

300.000  5000.000  1416.000 

5.31843737E-10-3. 04168815E-14 
4 . 904007 62E-02 -4 . 56067276E-05 
3 . 4318530  9E  +  01 

300.000  5000.000  1415.000 

6.59997414E-10-3. 78462012E-14 
5.26021422E-02-4 . 58866872E-05 
3 . 78560375E+01 

300.000  5000.000  1377.000 

8 . 58928005E-10-4 . 999844 92E-1 4 
2 . 94344471E-02-9. 79064273E-06 
1 . 60207425E+01 

300.000  5000.000  1398.000 

8 . 61 68014 6E-10-4 . 96213850E-14 
3. 93239957E-02-2 . 48355443E-05 
2 . 538818  91E  +  01 

300.000  5000.000  1399.000 

8 . 63922 637E-1 0-4 . 97707528E-14 
4 . 48313850E-02-3.21345956E-05 
2 . 245922  95E  +  01 

300.000  5000.000  1400.000 

8 . 65987714E-10-5 . 05377 999E-14 
5. 8130421 6E -02-4 . 457567 65E-05 
4 . 15138582E+01 

300.000  5000.000  1392.000 

8 . 45860293E-10-4 . 88215034E-14 
5. 87881110E-02-4 . 7 9917 613E-05 
3 . 02451397E+01 

300.000  5000.000  1418.000 

9. 1580181 9E-10-5 .25 63132 IE-14 
5. 94956538E-02-4 . 62659544E-05 
3 . 75273996E+01 

300.000  5000.000  1392.000 

1 . 09751772E-09-6. 32060952E-14 
4 . 0321203 6E -02-1 . 95983508E-05 
1 . 73968930E+01 

300.000  5000.000  1394.000 

1. 65880732E-09-9. 64168975E-14 
9. 60910608E-02-7 . 18927947E-05 
4.71199220E+01 

300.000  5000.000  1403.000 

1 . 19318832E-09-6. 904554 95E-14 
1 . 04258843E-01-9. 67113872E-05 
5. 92838342E+01 

300.000  5000.000  1397.000 

1 . 15824454E-09-6. 7443023 6E -14 
9.38131205E-02-9. 00665415E-05 
3 . 668584  99E+01 

300.000  5000.000  1400.000 

1 .20391708E-09-6. 98662618E-14 
1 . 0659 9512E-0 1-1 . 0583084 9E-0 4 
4 . 5181 9888E+01 

300.000  5000.000  1399.000 

1 .35877977E-09-7 . 9135 632 OE-14 
1 . 02109939E-01-8 . 86409603E-05 
5 . 74743892E+01 
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CVCZCVCJ  ZCVCECVC J_BK1 . C  8H  13N  OO  OG 

2 . 0  63027  83E  +  01  3 . 160  964 97E-02-1 . 0 69334 0  6E-05 
1 . 00684682E+04-7 . 7591054 9E+01-8 . 73 96664 4E-01 
1 . 67164961E-08-2 . 22051441E-12  1 . 77593242E+04 
CVCCCOJV  9/14/  9  THERMO  8H  130  1  OG 

2 . 4  94507  80E  +  01  3 . 0 660 6517E-02-1 . 05  90  657  3E-05 
-1. 695888 63E+03-1. 00312 800E+02-2. 8525 9996E+00 
2 . 97846865E-08-4 . 91192074E-12  7 . 7 60 60818E+03 
CJCVCCVC  CJCVCCVC_G3MP2C  9H  15N  OO  OG 

2 . 312  62384E  +  01  3 . 611 93184E-02-1 . 22113  663E-05 
1 .28522075E+04-8 . 62856416E+01  3 . 224 91856E+00 
1 . 0758167 6E -08-9. 72030195E-13  2 . 04256370E+04 
CVCCVCCOJ  9/14/  9  THERMO  9H  150  1  OG 

2 . 687512 63E+01  3 . 5982 6535E-02-1 . 22801111E-05 
-5 . 09466492E+03-1 . 0791645 6E +02-3 . 13932770E+00 
2 . 85244844E-08-4 . 3982 9458E-12  5 . 2 9230351E+03 
BZBCCLV  GA  C  12H  19  OCL  1G 

3 . 69272158E+01  4 . 4 6080 955E-02-1 . 52 93858 0E-05 
-1. 37 613553E+04-1. 5808 938 6E+02-2 . 4584 407 0E+00 
4 .24358334E-08-6. 9808012 8E-12-4 . 21163780E+02 
BZBBJ  9/  3/  9  THERMO  12H  19  0  OG 

3 . 34182780E+01  4 . 53424770E-02-1 . 55200912E-05 
7 . 80723881E+03-1 . 4188605 6E +02-5 . 06283930E+00 
3 . 94256080E-08-6 . 32734200E-12  2 . 09736691E+04 
BZBCCOJV  9/14/  9  THERMO  12H  190  1  OG 

3 . 582 6114 6E  +  01  4 . 55360593E-02-1 . 56116413E-05 
-3.26823346E+03-1 . 51 9334 64E+02-4 . 82372585E+00 
4 . 52844427E-08-7 . 54806140E-12  1 . 03936648E+04 
BBBCJV  9/  8/  9  THERMO  15H  23  0  OG 

4 . 21443837E+01  5 . 41028924E-02-1 . 83763736E-05 
9. 6465560 6E+03-1 . 847 98280E+02-5 . 41422303E+00 
4 . 58947967E-08-7 . 10170752E-12  2 . 59398952E+04 
BBBCOJV  9/14/  9  C  15H  230  1  OG 

4 . 69273361E+01  5 . 08221788E-02-1 . 68954273E-05 
-2 .24128532E+03-2 . 064 4757 1E+02-8 . 70352 630E+00 
6. 62988458E-08-1 . 11863416E-11  1 . 56309957E+04 
BBBJVIV  9/23/  9  C  16H  23  0  OG 

4 . 53398889E+01  5 . 427 97 942E-02-1 . 84165219E-05 
1 . 30903627E+04-2 . 02489154E+02-8 . 49422360E+00 
6. 70603758E-08-1 . 15989944E-11  3 . 05087 87 0E+04 
BBBVIV  9/21/  9  THERMO  16H  24  0  OG 

4 . 5530  982  9E  +  01  5 . 5  6991305E-02-1 . 87 92  7  585E-05 
1 . 0185456 6E +03-2 . 017 69691E+02-6 . 7477 9653E+00 
6 . 31731157E-08-1 . 0802055 9E -11  1 . 80887 113E+04 
R45MLSJ  GA  C  20H  32  OCL  1G 

5.81197  997E  +  01  7 . 7  02 98838E-02-2 . 5  907  4195E-05 
-1 . 08787019E+04-2 . 5627112 6E+02 -7 . 75715090E+00 
7 . 48361961E-08-1 . 23 66322 8E-11  1 . 08639815E+04 
R45MLPJ  GA  C  20H  32  OCL  1G 

5.692  92555E  +  01  7 . 7 6902385E-02-2 . 605  68  680E-05 
-8 . 31259134E+03-2 . 47362318E+02-5 . 1 60 97 675E+00 
7 . 11121480E-08-1 . 18094929E-11  1 . 22664116E+04 
R45M9J  9/14/  9  C  20H  31  0  OG 

5. 41954130E+01  7 . 60118 675E-02-2 . 57 1572 80E-05 
4 . 1428490 6E +03-2 . 423544 96E+02-8 . 4 9772505E+00 
6. 6164578 5E -08-1 . 06448625E-11  2 . 5305122 6E+04 
R45M13J  9/14/  9  C  20H  31  0  OG 

5 . 37073656E+01  7 . 64348507E-02-2 . 58637427E-05 
3. 1636353 6E+03-2. 3960596 6E+02-7 . 68140544E+00 
6.13997057E-08-9. 6232232 3E -12  2 . 41220366E+04 


300.000  5000.000  1399.000 

1 . 64774759E-09-9. 5092982 9E -14 
7 . 96074824E-02-5 . 0 935 98 94E-05 
3.86897  908E+01 

300.000  5000.000  1394.000 

1. 65509800E-09-9. 64 66438 6E-14 
9. 8122035 9E -02-7 . 45350560E-05 
4 . 81750535E+01 

300.000  5000.000  1393.000 

1 . 8807  031 9E-0  9-1 . 084  92  875E-13 
7 . 65942108E-02-4 . 17019077E-05 
2 . 282  93599E  +  01 

300.000  5000.000  1393.000 

1 . 903  93513E-0  9-1 . 103  65105E-13 
1 . 066072 04E-01-7 . 62737417E-05 
5 . 31236804E+01 

300.000  5000.000  1394.000 

2 . 3783205 6E -09-1 . 3814 630 9E-13 
1 . 40525666E-01-1 . 0 637 6675E-04 
5 . 21382  628E  +  01 

300.000  5000.000  1394.000 

2 . 41081354E-09-1 . 39924170E-13 
1 .37664846E-01-1 . 0152 114 0E-04 
6 . 39958354E+01 

300.000  5000.000  1396.000 

2 . 427 6272 6E -09-1 . 41003538E-13 
1 . 45587197E-01-1 . 118155 64E-04 
6 . 4  6655280E  +  01 

300.000  5000.000  1392.000 

2 . 8405545 6E -09-1 . 64330500E-13 
1 . 67056317E-01-1 . 214 66085E-04 
6 . 9902  61 62E  +  01 

300.000  5000.000  1395.000 

2 . 57613604E-09-1 . 47682896E-13 
1 . 9326340 6E -01-1. 57828780E-04 
8 . 75685214E+01 

300.000  5000.000  1397.000 

2 . 84443257E-09-1 . 6445304 9E-13 
1 . 925927 63E-01-1 . 57233014E-04 
8 . 2144  9680E  +  01 

300.000  5000.000  1395.000 

2 . 8 918453 9E-0 9-1 . 66770818E-13 
1 . 88529451E-01-1 . 50543355E-04 
7.517  96905E  +  01 

300.000  5000.000  1401.000 

3. 9770819 6E -09-2 . 2893791 9E -13 
2 . 41490030E-01-1 . 85182312E-04 
9 . 36459570E+01 

300.000  5000.000  1400.000 

3. 9921856 6E -09-2 . 2948235 9E -13 
2 . 32449762E-01-1 . 7 61 91351E-04 
8 . 2  630  9623E+01 

300.000  5000.000  1397.000 

3. 96336250E-09-2 . 287822 60E-13 
2 . 2 843754 7E-01-1 . 69268628E-04 
9 . 22049786E+01 

300.000  5000.000  1397.000 

3. 98670412E-09-2 . 3015154 3E -13 
2 . 23309632E-01-1 . 614182 68E-0 4 
8 . 88343536E+01 
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R45M  9/  3/  9  THERMO  20H  32  0  OG  300.000  5000.000  1398.000  91 

5 . 4250  657  8E  +  01  7 . 80424196E-02-2 . 63133838E-05  4 . 04 604518E-0 9-2 . 33170301E-13  2 

-1 .27596349E+04-2 . 41212500E+02-8 . 28123127E+00  2 . 30 6558 61E-01-1 . 70734012E-04  3 

6. 6953550 6E -08-1 . 0821814 9E-11  8 . 2 938 6157E+03  9 . 22860686E+01  4 

CL  42189CL  1  G  0300.00  5000.00  1000.00  1 

0 . 02920237E+02-0 . 03597 985E-02  0 . 012 942 94E-05-0 . 02162776E-09  0 . 0137 6517E-13  2 

0 . 01371338E+06  0 . 032 62 690E+02  0 . 02381577E+02  0 . 088 9107 9E-02  0 . 0407047 6E-05  3 

-0.02168943E-07  0 . Oil 60827E-10  0 . 01383999E+06  0 . 06021818E+02  4 

HCL  42189CL  1H  1  G  0300.00  5000.00  1000.00  1 

0 . 02755335E+02  0 . 01473581E-01-0 . 04 971254E-05  0 . 08108 658E-0 9-0 . 05072063E-13  2 

-0 . 01191806E+06  0 . 0 651511 6E  +  02  0 . 03338534E+02  0 . 012 68207E-01-0 . 03666917E-04  3 

0 . 04703992E-07-0 . 01836011E-10-0 . 01213151E+06  0 . 031 93555E+02  4 

CLOH  CLOH_G4 . out  CL  1H  IN  00  1G  300.000  5000.000  1425.000  01 

4 . 4381 961 9E+00  1 . 90159701E-03-5 . 96736743E-07  8 . 70883543E-11-4 . 82905289E-15  2 

-1 . 06618194E+04  2 . 38417393E+00  3 . 0 696232 9E+00  6 . 26906774E-03-5 . 88061081E-06  3 

2 . 91620508E-09-5 . 655104 63E-13-1 . 03069268E+04  9 . 312 95623E+00  4 

HOCLO  RE-FIT  0CL  10  2H  1  0G  300.000  4500.000  1000.00  0  1 

7 . 74227524E+00  1 . 14072335E-03-2 . 31095854E-07  3 . 98547731E-12  2 . 08972960E-15  2 

-1 . 45044 632E+02-1 .21025028E+01  2 . 33550334E+00  2 . 1821254 9E-02-3 . 3103904 6E-05  3 

2 . 5335102 9E -08-7 . 72997812E-12  1 . 0 92 69873E+03  1 . 44 913387E+01  4 

HOOCL  HOOCL_G4_5 0 0  0CL  10  2H  1  OG  300.000  5000.000  1000.00  01 

7 . 30200672D+00  1 . 611 92543D-03-2 . 7 90 65802D-07-8 . 62328785D-12  4 . 19722528D-15  2 

-3. 95942041 D+ 03-1. 06943293 D+ 01  2 . 04355669D+00  2 . 22 975388D-02-3 . 537 93621D-05  3 

2 . 91689375 D- 08-9. 50023046 D- 12-2 . 73718311D+03  1 . 51445456D+01  4 

H0CL02  RE-FIT  OCL  10  3H  1  OG  300.000  4500.000  1000.00  0  1 

9 . 81116772E+00  1 . 69027515E-03-3 . 3812 90 98E-07  6 . 20168803E-12  2 . 957 64 614E-15  2 

-4 . 73 63 6133E+03-2 .223110 96E+01  2 . 01147413E+00  3 . 2 0 615 61 9E-02-4 . 95894456E-05  3 

3 . 85937575E-08-1 . 190 4877 3E-1 1-2 . 97052759E+03  1 . 60160770E+01  4 

HCL04  RE-FIT  OCL  10  4H  1  OG  300.000  4500.000  1000.00  01 

1 . 15020361E+01  2. 5295130 9E -03-4. 6013431 6E -07-1. 21791145E-11  6 . 90475085E-15  2 

-4 . 85042 920E+03-3 . 10284100E+01  7 . 20484078E-01  4 . 04154 994E-02-5 . 86621500E-05  3 

4 . 58601157E-08-1 . 47 678085E-11-2 . 12432 935E+03  2 . 30658550E+01  4 

NH4CL04  RE-FIT  OCL  10  4N  1H  4G  300.000  4500.000  1000.00  01 

1 . 8317  8  654E  +  01  5 . 24443667E-03-8 . 23817 004E-07-5 . 82  0  90  62  6E-11  1 . 65  94  8001E-14  2 

-2 . 06336816E+04-6 . 58034286E+01  8 . 88326347E-01  6 . 68014288E-02-9 . 49331952E-05  3 

7 . 32077581E-08-2 . 32 674401E-11-1 . 627252 93E+0 4  2 . 14922791E+01  4 

AP_NH3  COMP_I RC  F_T  S_HC  L  1H  7N  20  4G  300.000  5000.000  1398.000  01 

1 . 82581064E+01  1 . 71061668E-02-5 . 61372775E-06  8 . 47432147E-10-4 . 82092675E-14  2 

-3. 3678863 6E +04-6. 80892417E+01  4 . 65325755E+00  5 . 32755888E-02-4 . 33406078E-05  3 

1.8 969507 7E-08-3. 3837 684 7E-12-2 . 93871425E+04  3 . 4303434 9E+00  4 

NOCL  CLNO_G4 . out  CL  IN  10  1  OG  300.000  5000.000  1378.000  01 

5 . 53178844E+00  1 . 3007 664 7E-03-4 . 61294870E-07  7 . 3372 9585E-11-4 . 32 941107E-15  2 

4 . 44818405E+03-8 . 99851231E-01  4 . 13186473E+00  4 . 51 65 67 94E-03-3 . 33281748E-06  3 

1 .26992097E-09-2 . 01164271E-13  4 . 9456512 6E+03  6 . 64809180E+00  4 

0N0CL03  RE-FIT  OCL  IN  10  5  OG  300.000  4500.000  1000.00  01 

1. 67813644E+01  1 . 11842 831E-03-3 . 38 9212 68E-07  3 . 39122758E-11-1 . 90400672E-16  2 

1 . 27 91710 9E +04-5 . 48439751E+01  3 . 101 92537E+00  4 . 2 990177 9E-02-5 . 48648932E-05  3 

3 . 81883822E-08-1 . 18209947E-11  1 . 6536142 6E+04  1 . 52785330E+01  4 

CLO  RE-FIT  OCL  10  1  0  OG  300.000  4500.000  1000.00  0  1 

4 . 23923159D+00  1 . 9344 92 61D-04-4 . 25335251D-08  7 . 095512 87D-13  4 . 32 92 634 8D-1 6  2 

1 . 07 95014 6D+04  2 . 1514  7  68  6D+00  2 . 92  077  8  99D+00  4 . 0  90 4  54  4 3D-03-4 . 4 92  57  32  3D-0 6  3 

2 . 48358267D-09-6. 10953779D-13  1 . 114 98711D+04  8 . 9122 61 96D+00  4 

CL02  CL02_G4 . out  CL  1H  ON  OO  2G  300.000  5000.000  1408.000  01 

6 . 13221573E+00  8 . 17125271E-04-3 . 00956329E-07  4 . 90756142E-11-2 . 94593271E-15  2 

9. 74125871E+03-4 . 93151693E+00  2 . 78920691E+00  1 . 03639661E-02-1 . 07452184E-05  3 

5.19247877E-09-9.57201985E-13  1 . 07053068E+04  1 . 23677689E+01  4 

CLOO  CLOO_G4 . out  CL  IN  OO  2  OG  300.000  5000.000  1682.000  01 

5. 9856702 6E+00  9 . 7 630594 9E-04-3 . 67418941E-07  6 . 098237 60E-11-3 . 7 105527 8E-15  2 

1 . 02856449E+04-8 . 31833092E-01  5 . 06381486E+00  2 . 94867142E-03-1 . 81698983E-06  3 

4 . 80823618E-10-4 . 025577 47E-14  1 . 06037780E+04  4 . 15563337E+00  4 
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CL03  CL03_G4 . out  CL  1H  ON  OO  3G  300.000  5000.000  1407.000  01 

8 . 55671903E+00  1 . 34 914053E-03-4 . 94 913928E-07  8 . 0508857 6E-11-4 . 82528654E-15  2 

1. 9148679 9E +04-1. 77737893E+01  1 . 91755743E+00  2 . 152392 68E-02-2 . 39015577E-05  3 

1 .22170196E-08-2 . 35 60330 9E-12  2 . 0 954 6638E+04  1 . 61728591E+01  4 

CL2  42189CL  2  G  0300.00  5000.00  1000.00  1 

0 . 04274587E+02  0 . 03717337E-02-0 . 01893490E-05  0 . 05337 4 65E-0 9-0 . 05057 602E-13  2 

-0 . 01331149E+05  0 . 02256947E+02  0 . 03439587E+02  0 . 02870774E-01-0 . 02385871E-04  3 

0 . 02892918E-08  0 . 02 915057E-11-0 . 01131787E+05  0 . 0 64 7 135 9E+02  4 

H  120186H  1  G  0300.00  5000.00  1000.00  1 

0 . 02500000E+02  0 . 00000000E+00  0 . 00000000E+00  0 . 00000000E+00  0 . 00000000E+00  2 

0 . 02547163E+06-0 . 04 60117 6E+01  0 . 02500000E+02  0 . 00000000E+00  0 . 00000000E+00  3 

0.00000000E+00  0.00000000E+00  0 . 02547163E+06-0 . 0460117 6E +01  4 

HNO  WRA032 4 98  H  IN  10  1  0G  300.00  5000.00  1000.000  01 

3 . 16554762E+00  3 . 00005132E-03-3 . 94350282E-07-3 . 857874 91E-11  7 . 08091931E-15  2 

1 . 18052184E+04  7 . 647 64 695E+00  4 . 53525882E+00-5 . 6854 6910E-03  1 . 851 9997 6E-05  3 

-1 . 71883674E-08  5 . 558330 90E-12  1 . 16506820E+04  1 . 74314734E+00  1 . 28824477E+04  4 

HONO  31787H  IN  10  2  G  0300.00  5000.00  1000.00  1 

0 . 054  8  68  93E  +  02  0 . 042180  65E-01-0 . 0164  914  3E-04  0 . 02 97187  7E-08-0 . 02  02114  8E-12  2 

-0 . 01126865E+06-0 . 02997002E+02  0 . 022 90413E+02  0 . 01409922E+00-0 . 01367872E-03  3 

0 . 07  4 987  80E-07-0 . 0187 6905E-10-0 . 010431 95E+0 6  0 . 01328077E+03  4 

NNH  120186N  2H  1  G  0250.00  4000.00  1000.00  1 

0 . 04415342E+02  0 . 01 614388E-01-0 . 01632894E-05-0 . 08559846E-09  0 . 01614791E-12  2 

0 . 027  8802  9E  +  0  6  0 . 09042888E+01  0 . 03501344E+02  0 . 02053587E-01  0 . 07170410E-05  3 

0 . 04921348E-08-0 . 09671170E-11  0 . 02 83334 7E+0 6  0 . 06391837E+02  4 

HNJN02  HN JN02_G4 . chk  H  IN  20  2  0G  300.000  5000.000  1393.000  11 

8 . 5882  9360E  +  00  3 . 557 83 685E-03-1 . 37988061E-06  2 . 341402 97E-10-1 . 4205744 6E-14  2 

2 . 48692567E+04-1 . 7281030 9E +01  2 . 38719149E+00  2 . 0 6880 667E-02-1 . 97118815E-05  3 

9. 16137562E-09-1 . 66604445E-12  2 . 67339286E+04  1 . 504518 67E+01  4 

HNOJN02  HNO JN02_B3d_B3H  IN  20  3  0G  300.000  5000.000  1392.000  01 

1 . 03253946E+01  4 . 98 6752 64E-03-1 . 80524376E-06  2 . 92271754E-10-1 . 73600960E-14  2 

1 . 02125162E+04-2 . 57969338E+01  2 . 7 6417845E+00  2 . 48231368E-02-2 . 208427 98E-05  3 

9. 8134383 9E -09-1 . 73206394E-12  1 . 26074164E+04  1 . 40188663E+01  4 

OH  Sand/RuscOl  0  1H  1  G  0300.00  5000.00  1000.00  1 

0 . 02882730E+02  0 . 01013974E-01-0 . 0227 6877E-05  0 . 02174684E-09-0 . 05126305E-14  2 

0 . 03672089E+05  0 . 055 957 12E+02  0 . 036372 66E+02  0 . 01850 910E-02-0 . 01 67 61 65E-04  3 

0 . 02387203E-07-0 . 08431442E-11  0 . 033 91 983E+05  0 . 013588 60E+02  4 

H02  McB 93/ WRA0 0  H  10  2  0  OG  300.00  5000.00  1000.000  01 

4 . 17228728E+00  1 . 88117 64 7E-03-3 . 4 6277408E-07  1 . 94657853E-11  1 . 7 62542 94E-1 6  2 

2 . 134177 95E  +  02  2 . 957 6774 6E+00  4 . 3017 9801E+00-4 . 74 912051E-03  2 . 11582891E-05  3 

-2.42763894E-08  9 . 2 9225124E-12  4 . 4 64 1553 9E+02  3 . 71666245E+00  1 . 66125750E+03  4 

H2  121286H  2  G  0300.00  5000.00  1000.00  1 

0 . 02 991423E  +  02  0 . 07 000 64 4E-02-0 . 05 63382 9E-0 6-0 . 0  92  3157 8E-10  0 . 01582752E-13  2 

-0 . 0 835034 0E+0 4-0 . 0 13551 10E+02  0 . 03298124E+02  0 . 0824 9442E-02-0 . 08143015E-05  3 

-0 . 09475434E-09  0 . 04 134 872E-11-0 . 0101252 1E+05-0 . 032 940 94E+02  4 

NH2  McB93  N  1H  2  0  OG  300.00  5000.00  1000.000  01 

2.847 68992E  +  00  3 . 142  80035E-03-8 . 98  641458E-07  1 . 303182 84E-10-7 . 4  8812 92 6E-15  2 

2 . 1823904 9E  +  04  6 . 47165433E  +  00  4 . 20556857E  +  00-2 . 13561363E-03  7 . 2 6851301E-0 6  3 

-5. 9306987 6E -09  1 . 80690978E-12  2 . 15352231E+04-1 . 4 6662770E-01  2 . 27475415E+04  4 

NH20  M/JB86  N  1H  20  1  OG  300.000  5000.000  1398.000  01 

4 . 2 6222 939E  +  00  4 . 60071183E-03-1 . 52686779E-06  2 . 32081624E-10-1 . 32607907E-14  2 

6 . 2  6937 941E  +  03  1 . 8 952  3882E  +  00  2 . 62132  814E  +  00  8 . 055 942 93E-03-4 . 341997  52E-0  6  3 

1 . 31067 689E-09-1 . 7 941316 9E -13  6 . 89825870E  +  03  1 . 087 68221E  +  01  4 

N2H2  RB91/ Gurvich  N  2H  2  0  OG  300.00  5000.00  1000.00  01 

1 . 43773356E+00  8 . 72166734E-03-2 . 99323466E-06  4 . 5002 9732E-10-2 . 50422 693E-14  2 

2 . 22727357E+04  1 . 574532 60E+01  4 . 74387065E+00-9 . 15242446E-03  3 . 3534 98 97E-05  3 

-3.22086920E-08  1 . 06734098E-11  2 . 17525597E+04  7 . 28947240E-01  2 . 29978160E+04  4 

H2NN02  RE-FIT  OH  2N  20  2  OG  300.000  4500.000  1000.00  01 

9 . 45441628E+00  3 . 80055304E-03-5 . 10139330E-07-7 . 9707 9347E-11  1 . 6241534 9E-14  2 

-3. 15945020E+03-2. 5396232 6E +01  4 . 51412708E-01  3 . 07305455E-02-3 . 94054223E-05  3 

3. 2932902 6E -08-1. 2028075 6E -11-5. 00637573E+02  2 . 13188362E+01  4 
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H20  20387H  20  1  G  0300.00  5000.00  1000.00  1 

0 . 02 67214 6E  +  02  0 . 030562 93E-01-0 . 087302 60E-05  0 . 01200996E-08-0 . 06391618E-13  2 

-0 . 02989921E+06  0 . 06862817E+02  0 . 0338 6842E+02  0 . 03474 982E-01-0 . 06354696E-04  3 

0 . 06968581E-07-0 . 02506588E-10-0 . 03020811E+06  0 . 02590233E+02  4 

NH3  121386N  1H  3  G  0300.00  5000.00  1000.00  1 

0 . 02461904E+02  0 . 0605 9166E-01-0 . 02004 97 7E-04  0 . 03136003E-08-0 . 01938317E-12  2 

-0 . 06493270E+05  0 . 07472097E+02  0 . 02204352E+02  0 . 0101147 6E+00-0 . 014 652 65E-03  3 

0. 014 47235E-0 6-0. 0532 850 9E-10-0.0 65254 88E+05  0 . 08127 138E+02  4 

N2H3  DB00/RB91  N  2H  3  0  0G  300.000  5000.000  1385.000  11 

5 . 56523220E+00  5 . 92789098E-03-2 . 05378972E-06  3 . 21563158E-10-1 . 87 648230E-14  2 

2 . 40692822E+04-7 . 11150205E+00  1 . 68398813E+00  1 . 4 64554 61E-02-9 . 84764677E-06  3 

3. 66967922E-09-5. 99668094E-13  2 . 54 978433E+04  1 . 39338819E+01  2 . 65708384E+04  4 

0JHNH3  0 JHNH3_BK1_G4  CL  OH  4N  10  1G  300.000  5000.000  1417.000  01 

6.1517 6103E+00  7 . 14174114E-03-2 . 20334210E-06  3 . 17596288E-10-1 . 74502751E-14  2 

-7 . 40333472E+03-8 . 02764592E+00  4 . 0843573 6E+00  1 . 2 957 62 97E-02-8 . 64359256E-06  3 

3. 57714495E-09-6. 4162562 6E -13-6. 77 69209 6E +03  2 . 73636599E+00  4 

N2H4  B0Z_JPCA1 1 4_62H  4N  2  0  OG  300.000  5000.000  1390.000  11 

5 . 02599050E+00  8 . 65302 984E-03-2 . 93055474E-06  4 . 51741564E-10-2 . 60716879E-14  2 

8. 97814993E+03-3. 92585332E+00  1 . 56439782E+00  1 . 55016475E-02-8 . 06079743E-06  3 

2 . 24077 977E-09-2 . 80712108E-13  1 . 03609178E+04  1 . 51984751E+01  4 

DINH3  DINH3_BK1_1R0TCL  OH  6N  20  OG  300.000  5000.000  1429.000  11 

6 . 21669185E+00  1 . 12895910E-02-3 . 52568970E-06  5 . 12341681E-10-2 . 83064658E-14  2 

-1 . 4234484 6E +04-1 .29402375E+00  6 . 7 942 90 66E+00  8 . 88257 950E-03-5 . 13900006E-07  3 

-1 . 00249924E-09  2 . 4322592 6E-13-1 . 427 6217 OE+O 4-3 . 93519667E+00  4 

NO  MCB93  N  10  1  0  OG  300.00  5000.00  1000.000  01 

3 . 2 6071234E  +  00  1 . 1 9101135E-03-4 . 2 9122 64 6E-07  6 . 944814 63E-11-4 . 032 95681E-15  2 

9. 92143132E+03  6 . 36900518E+00  4 . 21859896E+00-4 . 63988124E-03  1 . 1044304 9E-05  3 

-9.34055507E-09  2 . 80554874E-12  9 . 8450 9964E+03  2 . 28061001E+00  1 . 09770882E+04  4 

N02  121286N  10  2  G  0300.00  5000.00  1000.00  1 

0 . 04 682859E  +  02  0 . 02  4 62  42 9E-0 1-0 . 0104225 9E-04  0 . 01 97 6902E-08-0 . 013  91717E-12  2 

0.022 612 92E+05  0 . 09885985E+01  0 . 02670600E+02  0 . 07838501E-01-0 . 08063865E-04  3 

0 . 06161715E-07-0 . 02320150E-10  0 . 02896291E+05  0 . 01161207E+03  4 

N2  121286N  2  G  0300.00  5000.00  1000.00  1 

0 . 02926640E+02  0 . 01487 977E-01-0 . 056847 61E-05  0 . 01009704E-08-0 . 06753351E-13  2 

-0 . 09227 977E+04  0 . 05 98052 8E+02  0 . 032 98 67 7E+02  0 . 01408240E-01-0 . 03963222E-04  3 

0. 05641515E-07-0.02444855E-10-0. 01020 900E+05  0 . 03 950372E+02  4 

0  1201860  1  G  0300.00  5000.00  1000.00  1 

0 . 02542060E+02-0 . 02755062E-03-0 . 03102803E-07  0 . 045510 67E-1 0-0 . 04368052E-14  2 

0 . 02 923080E  +  0 6  0 . 04 920308E  +  02  0 . 02 94 642 9E  +  02-0 . 01 6381 66E-01  0 . 02421032E-04  3 

-0.01602843E-07  0 . 038 90 696E-11  0 . 02914764E+06  0 . 02 963995E+02  4 

02  1213860  2  G  0300.00  5000.00  1000.00  1 

0 . 03697578E+02  0 . 0 61351 97E-02-0 . 01258842E-05  0 . 01775281E-09-0 . 01136435E-13  2 

-0 . 01233930E+05  0 . 0318 91 66E+02  0 . 032 12 93 6E+02  0 . 01127486E-01-0 . 05756150E-05  3 

0. 0131387 7E-07-0. 087 68554E-11-0. 010052 4 9E+05  0 . 06034 738E+02  4 

END 
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Intentionally  lelt  blank. 
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Appendix  C.  The  Molecular  Structures  of  Species  in 

Mechanism  R2 
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Because  the  CHEMKIN  II  framework  limits  the  labels  for  species  to  16  or  fewer 
characters  and  many  of  the  species  are  radicals  or  “unusual”,  in  most  cases  it  is  not 
possible  to  assign  a  label  (such  as  a  SMILES  [simplified  molecular-input  line-entry 
system]  string)  that  could  be  unambiguously  translated  into  a  molecular  structure. 
Correspondences  between  the  labels  and  the  molecular  structures  for  which  the 
stoichiometry  and  (single)  mandate  that  hydrogen  atoms  only  form  single  bonds  are 
not  sufficient  to  yield  an  unambiguous  interpretation  are  provided  in  Table  C-l. 


Table  C-l  The  molecular  structures  of  species  in  mechanism  R2a 


Label 

Stoichiometry 

Structure 

CHO 

CHO 

CH20 

CH20 

y  % 

C02 

CO2 

§-i-« 

C2H2 

C2H2 

C2H3 

2c2h3 

C*CCL*0 

C3C1H30 

VCDJO 

2c3h3o 

Jyr* 

J 

VCDO 

C3H40 

XX* 

J 

CTCV 

C4H4 

r 

CDCCDCJ 

2c4h5 

j 
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CDCCDC 

c4h6 

CCDOV 

c4h6o 

CDCCJC 

2c4h7 

VCOJC 

2c4h7o 

CJVV 

2c5h7 

WCOJ 

2c5h7o 

vVr 

CVV 

c5h8 

jfyls JL* 

J  J 

CJCVC 

2c5h9 

CVCCVCL 

CsCIHn 

VY6DE13J 

2c8h9 

& 
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VVVJV 

2c8h9 

VVVIV 

CgHio 

VrV- 

VY6DE13 

CsHio 

CVCZCVCJ 

2C8Hl3 

* 

CVCCCOJV 

2C8Hi30 

CJCVCCVC 

2C9Hi5 

CVCCVCCOJ 

2C9Hi50 

BZBCCLV 

Ci2C1Hi9 

BZBBJ 

2Ci2Hi9 
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BZBCCOJV 

2Ci2Hi90 

BBBCJV 

2Ci5H23 

■H 

BBBCOJV 

2Ci5H230 

BBBJVIV 

2Ci6H23 

V* 

T 

BBBVIV 

Ci6H24 

A 
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R45MLSJ 

2C2oC1H32 

R45MLPJ 

2C2oC1H32 

R45M9J 

2C2oH3i 

R45M13J 

2C20H3i 

¥ 

R45M 

C2oH32 

CLOH 

C1HO 

HOCLO 

CIHO2 

•  ^ 

HOOCL 

CIHO2 
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H0CL02 

C1H03 

HCL04 

CIHO4 

v* 

NH4CL04 

CIH4NO4 

APNH3 

CIH7N2O4 

NOCL 

C1NO 

•  •* 

0N0CL03 

CINO5 

CL02 

2C102 

CLOO 

2cio2 

•  % 

CL03 

2cio3 

•r 

HNO 

HNO 

HONO 

hno2 

HNJN02 

2hn2o2 

>♦. 

HN0JN02 

2hn2o3 

H 
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NH20 

2H2NO 

N2H2 

H2N2 

J 

H2NN02 

H2N2O2 

K 

N2H3 

2h3n2 

W 

J 

OJHNH3 

2h4no 

4  - 

N2H4 

H4N2 

'VC 

DINH3 

h6n2 

4 

acarbon  (C)  atoms  =  gray;  hydrogen  (H)  atoms  =  white;  nitrogen  (N)  atoms  =  blue;  oxygen  (O) 
atoms  =  red;  and  chlorine  (Cl)  atoms  =  green. 
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Appendix  D.  Parameters  Employed  For  Computing  Transport 
Property  Estimates  for  Species  in  Mechanism  R2 
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Table  D-l  provides  values  for  the  coefficients  employed  to  estimate  the  transport  properties  of  pure  species  in  mechanism  R2.  For  each 
species/molecule,  there  is  a  set  that  includes  1)  an  index  indicating  whether  it  is  monatomic,  linear,  or  nonlinear,  2)  its  Lennard- Jones 
potential  well  depth  (s/fo),  3)  its  Lennard-Jones  collision  diameter  (cr),  4)  its  dipole  moment  (//),  5)  its  polarizability  (a),  and  6)  its 
rotational  relaxation  collision  number  (Zrot). 

Table  D-l  Transport  parameters  for  species  in  mechanism  R2 


Label 

8/ks 

G 

JU 

a 

Z rot 

CHO 

2 

71.4 

3.798 

0 

0 

1 

CH20 

2 

498 

3.590 

0 

0 

2 

CH3 

1 

144 

3.800 

0 

0 

0 

CO 

1 

98.1 

3.650 

0 

1.95 

1.8 

C02 

1 

244 

3.763 

0 

2.65 

2.1 

C2H2 

1 

209 

4.100 

0 

0 

2.5 

C2H3 

2 

209 

4.100 

0 

0 

1 

C*CCL*0 

2 

445 

4.890 

0 

0 

1 

VCDJO 

2 

425 

4.970 

0 

0 

1 

VCDO 

2 

425 

4.970 

0 

0 

1 

CTCV 

2 

357 

5.180 

0 

0 

1 

CDCCDCJ 

2 

357 

5.176 

0 

0 

1 

CDCCDC 

2 

357 

5.176 

0 

0 

1 

CCDOV 

2 

428 

5.320 

0 

0 

1 

CDCCJC 

2 

357 

5.176 

0 

0 

1 

VCOJC 

2 

428 

5.320 

0 

0 

1 

CJVV 

2 

348 

5.650 

0 

0 

1 

VVCOJ 

2 

473 

5.630 

0 

0 

1 

CVV 

2 

362 

5.470 

0 

0 

1 

CJCVC 

2 

362 

5.470 

0 

0 

1 

Remarks/Reference _ 

CHO:  Chemkin  (same  as  HCO) 

CH20:  Chemkin 
CH3:  Chemkin 
CO:  Chemkin 
C02:  Chemkin 
C2H2:  Chemkin 
C2H3:  Chemkin 

C3C1H30:  Assume  similar  to  CH3CH2N02  (Tb=387K;  Vb=71  cm3/mol) 
C3H30:  Assume  similar  to  1 -propanol  (Tb=370K;Vb=75  cm3 /mol) 
C3H40:  Assume  similar  to  1 -propanol  (Tb=370K;Vb=75  cm3 /mol) 

C4H4:  (=  C4H4  in  Chemkin) 

C4H5:  (=  C4H6  in  Chemkin) 

C4H6:  Chemkin 

C4H60:  Assume  similar  to  2-butanol  (Tb=372K;Vb=92  cm3/mol) 

C4H7  (=C4H6  in  Chemkin) 

C4H70:  Assume  similar  to  2-butanol  (Tb=372K;Vb=92  cm3/mol) 

C5H7:  Assume  similar  to  1-pentene  (Tb=303K;Vb=110  cm3  /mol) 

C5H70:  Assume  similar  to  1-pentanol  (Tb=41  lK;Vb=109  cm3/mol) 
C5H8:  1,3-pentadiene  (Tb=315K;Vb=99.7  cm3/mol) 

C5H9:  Assume  similar  to  1,3-pentadiene  (Tb=315K;Vb=99.7  cm3/mol) 
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Table  D-l  Transport  parameters  for  species  in  mechanism  R2  (continued) 


Label 

s/Ub 

G 

V 

a 

Z rot 

Remarks/Reference 

CVCCVCL 

2 

524 

6.540 

0 

0 

1 

C8C1H13:  1-chlorooctane  (Tb=456K;Vb=170  cm3/mol) 

VY6DE13J 

2 

470 

5.860 

0 

0 

1 

C8H9:  Assume  similar  to  ethylbenzene(Tb=409K;Vb=122  cm3  /mol) 

VVVIV 

2 

451 

6.440 

0 

0 

1 

C8H10:  Assume  similar  to  3-methylheptane  (Tb=392K;Vb=162  cm3/mol) 

VY6DE13 

2 

470 

5.860 

0 

0 

1 

C8H10:  Assume  similar  to  ethylbenzene(Tb=409K;Vb=122  cm3/mol) 

CVCZCVCJ 

2 

458 

6.440 

0 

0 

1 

C8H13:  Assume  similar  to  octane  (Tb=399K;Vb=162  cm3/mol) 

CVCCCOJV 

2 

520 

6.390 

0 

0 

1 

C8H130:  Assume  similar  to  2-octanol  (Tb=452K;Vb=159  cm3/mol) 

CJCVCCVC 

2 

488 

6.650 

0 

0 

1 

C9H15:  Assume  similar  to  nonane  (Tb=424K;Vb=179  cm3/mol) 

CVCCVCCOJ 

2 

560 

6.590 

0 

0 

1 

C9H150:  Assume  similar  to  1-nonanol  (Tb=487K;Vb=174  cm3/mol) 

BZBCCLV 

2 

601 

7.170 

0 

0 

1 

C12C1H19:  Assume  similar  to  2-dodecanol  (Tb=523K;Vb=225  cm3 /mol) 

BZBBJ 

2 

564 

7.200 

0 

0 

1 

C12H19:  Assume  similar  to  dodecane  (Tb=490K;Vb=227  cm3 /mol) 

BZBCCOJV 

2 

601 

7.170 

0 

0 

1 

C12H190:  Assume  similar  to  2-dodecanol  (Tb=523K;Vb=225  cm3/mol) 

BBBCJV 

2 

624 

7.530 

0 

0 

1 

C15H23:  Assume  similar  to  pentadecane  (Tb=543K;Vb=260  cm3  /mol) 

BBBCOJV 

2 

658 

7.680 

0 

0 

1 

C15H230:  Assume  similar  to  1 -pentadecanol  (Tb=572K;Vb=276  cm3/mol) 

BBBJVIV 

2 

644 

7.850 

0 

0 

1 

C16H23:  Assume  similar  to  cetane  (Tb=560K;Vb=294  cm3 /mol) 

BBBVIV 

2 

644 

7.850 

0 

0 

1 

C16H24:  Assume  similar  to  cetane  (Tb=560K;Vb=294  cm3/mol) 

R45MLSJ 

2 

742 

8.140 

0 

0 

1 

C20C1H32:  Assume  similar  to  icosanol  (Tb=645K;Vb=355  cm3/mol) 

R45MLPJ 

2 

742 

8.140 

0 

0 

1 

C20C1H32:  Assume  similar  to  icosanol  (Tb=645K;Vb=355  cm3 /mol) 

R45M9J 

2 

708 

8.380 

0 

0 

1 

C20H31:  Assume  similar  to  icosane  (Tb=616K;Vb=358  cm3/mol) 

R45M13J 

2 

708 

8.380 

0 

0 

1 

C20H31:  Assume  similar  to  icosane  (Tb=616K;Vb=358  cm3 /mol) 

R45M 

2 

708 

8.380 

0 

0 

1 

C20H32:  Assume  similar  to  icosane  (Tb=616K;Vb=358  cm3/mol) 

CL 

0 

130.8 

3.613 

0 

0 

1 

Cl:  Chemkin  2002 

HCL 

1 

344.7 

3.339 

0 

0 

1 

C1H:  Chemkin  2002 

CLOH 

2 

430 

3.690 

1.844 

0 

1 

C1HO:  Tb(est)=374K;  Tb(H20)=373K  (McQuaid  and  Chen  2014) 

HOCLO 

2 

469 

3.920 

1.844 

0 

1 

C1H02:  Tb(est)=408K  (McQuaid  and  Chen  2014) 

HOOCL 

2 

469 

3.920 

1.844 

0 

1 

C1H02:  Assume  same  as  HOCIO  (McQuaid  and  Chen  2014) 

H0CL02 

2 

508 

4.120 

1.844 

0 

1 

C1H03:  Tb(est)=442K  (McQuaid  and  Chen  2014) 
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Table  D-l  Transport  parameters  for  species  in  mechanism  R2  (continued) 


Label 

s/Ub 

G 

V 

a 

Z rot 

Remarks/Reference 

HCL04 

2 

547 

4.310 

1.844 

0 

1 

C1H04:  Tb(meas)=476K1 

NH4CL04 

2 

547 

4.770 

0 

0 

1 

C1H4N04:  Tb(est)=476K=Tb(HCL04,  est)  (McQuaid  and  Chen  2014) 

APNH3 

2 

547 

4.770 

0 

0 

1 

C1H7N204:  Assume  similar  to  NH4C104 

NOCL 

2 

450 

3.880 

0 

0 

1 

C1NO:  Svelha  1995;  delta=0.50 

ONOCL03 

2 

547 

4.770 

0 

0 

1 

C1N05:  Assume  similar  to  NH4C104 

CLO 

1 

288 

3.690 

0 

0 

1 

CIO:  Tb(est)=250K;  Tb(SO)=262  K;Tb(C12)=239  K;  Vb=Vb(CLOH) 

CL02 

2 

327 

3.920 

0 

0 

1 

C102:  Tb(meas)=284K;Vb=Vb(HOCLO) 

CLOO 

2 

308 

3.920 

0 

0 

1 

C102:  Tb(est)=268K  =  Tb(NOCL);  Vb=V(HOCLO) 

CL03 

2 

366 

4.120 

0 

0 

1 

C103:  Tb(est)=318  K;  Tb(S03)=318  K;  Vb=V(HOCL02) 

CL2 

1 

300 

4.235 

0 

0 

1 

C12:  Svelha  1995 

H 

0 

145 

2.050 

0 

0 

0 

H:  Chemkin 

HNO 

2 

116.7 

3.492 

0 

0 

1 

HNO:  Chemkin 

HONO 

2 

572.4 

3.600 

1.844 

0 

3 

HN02:  from  Penn  State  via  Miller  (combines  parameters  of  H2O  and  NO2) 

NNH 

2 

71.4 

3.798 

0 

0 

1 

HN2:  Chemkin 

HNJN02 

2 

409 

4.090 

0 

0 

1 

HN202:  Assume  similar  to  HN03  (Tb=356K;  Vb=42  cm3/mol) 

HNOJN02 

2 

389 

4.720 

0 

0 

1 

HN203:  Assume  similar  to  CH30N02  (Tb=338K;  Vb=64  cm3/mol) 

OH 

1 

80 

2.750 

0 

0 

0 

HO:  Chemkin 

H02 

2 

107.4 

3.458 

0 

0 

1 

H02:  Chemkin 

H2 

1 

38 

2.920 

0 

0.79 

280 

H2:  Chemkin 

NH2 

2 

80 

2.650 

0 

2.26 

4 

H2N:  Chemkin 

NH20 

2 

116.7 

3.492 

0 

0 

1 

H2NO:  Assume  similar  to  HNO 

N2H2 

2 

71.4 

3.798 

0 

0 

1 

H2N2:  Chemkin 

H2NN02 

2 

430 

3.810 

0 

0 

1 

H2N202:  Nitramide.  Tb(est)=374K;  Vb(Klapotke)=34  cm3/mol) 

H20 

2 

572.4 

2.605 

1.844 

0 

4 

H20:  Chemkin 

NH3 

2 

481 

2.920 

1.47 

0 

10 

H3N:  Chemkin 

N2H3 

2 

200 

3.900 

0 

0 

1 

H3N2:  Chemkin 

OJHNH3 

2 

381 

3.550 

0 

0 

1 

H4NO:  Assume  similar  to  NH20H  (Tb=331K;  Vb=27.6  cm3/mol) 

^cQuaid  MJ,  Chen  CC.  Modeling  the  deflagration  of  ammonium  perchlorate  at  pressures  from  300  to  30000  psia.  Part  II:  considerations  besides  the  gas-phase,  finite-rate  chemical 
kinetics  mechanism.  Presented  at  the  46th  JANNAF  Combustion  Subcommittee  Meeting;  2014  Dec  8-11;  Albuquerque,  NM. 
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Table  D-l  Transport  parameters  for  species  in  mechanism  R2  (continued) 


Label 

s/Ub 

G 

V 

a 

Z rot 

Remarks/Reference 

N2H4 

2 

205 

4.230 

0 

4.26 

1.5 

H4N2:  Chemkin 

DINH3 

2 

276 

4.350 

0 

0 

1 

H6N2:  Assume  Tb=Tb(NH3)=240K;  Vb=2*Vb(NH3)=50  cm3/mol 

NO 

1 

97.53 

3.621 

0 

1.76 

4 

NO:  Chemkin 

N02 

2 

200 

3.500 

0 

0 

1 

N02:  Chemkin 

N2 

1 

97.53 

3.621 

0 

1.76 

4 

N2:  Chemkin 

0 

0 

80 

2.750 

0 

0 

0 

O:  Chemkin 

02 

1 

107.4 

3.458 

0 

1.6 

3.8 

02:  Chemkin 

Except  for  the  index  indicating  whether  a  molecule  is  monatomic,  linear,  or 
nonlinear,  for  all  but  common  combustion  products,  values  for  these  parameters  are 
not  obvious  or  obtainable  from  open  literature  sources  and  have  to  be  estimated. 
While  this  has  the  potential  to  be  an  arduous  task,2  we  found  negligible  differences 
in  burning  rate  estimates  for  pure  ammonium  perchlorate  produced  when  a  set  for 
a  relatively  small  and  nonpolar  molecule  such  as  CH3CHO  was  employed  for  all 
unknowns  and  when  the  sets  included  first-order  approximations  for  e/fo  and  o  and 
a  value  for  n  =  1.844  (instead  of  zero)  if  the  molecule  had  a  hydroxyl  group 
(McQuaid  and  Chen  2014).  Concluding  that  obtaining  first-order  approximations 
for  dks  and  a  was  adequate  for  our  purposes,  we  employed  that  approach  to 
estimate  values  for  species  in  mechanism  R2. 

For  species  for  which  established  estimates  for  s/Icb,  o,  h,  a,  and  Zrot  could  not  be 
found,  we  followed  past  practice:  setting  a  =  0  and  Zrot  =  1.0.  If  the  molecule  had 
a  hydroxyl  group,  n  was  set  equal  to  1.844  (i.e.,  the  //  for  H2O).  Otherwise  it  was 
set  equal  to  0. 

Estimates  for  e/Icb  and  a  were  obtained  via 

s/kB=  1.15*  Tb , 

where  Tb  is  the  species’  normal  boiling  point  and 

a=\.lS*Vbm, 


where  V/,  is  the  species’  molecular  volume  at  the  normal  boiling  point  (Welty  et  al. 
1976). 3  Values  for  Tb  and  V/,  were  obtained  via  an  internet  search.  Wikipedia  was 
typically  consulted  first.  If  it  did  not  yield  a  suitable  result,  the  Sigma-Aldrich 
catalog  and/or  www.nbchem.ncbi.nlm.nih.gov  were  searched.  None  of  the  Tb  and 
Vb  values  that  were  found  had  formal  attribution.  Our  only  criteria  for  employing  a 
value  was  that  it  conform  to  our  expectations  for  the  relative  values  of  species  based 
on  molecular  size  and  the  presence  of  certain  functional  groups.  Moreover,  in  the 
majority  of  cases  (including  most  radicals),  Tb  and  Vb  for  the  molecule/species  of 
specific  interest  could  not  be  (nor  could  they  be  expected  to  be)  found.  In  these 
cases,  values  for  a  species  of  similar  elemental  composition/size  and  polarity  were 
sought  and  used. 


9 

Troe,  J.  Theory  of  thermal  unimolecular  reactions  at  low  pressures:  II.  strong  collision  rate  constants: 
applications.  Journal  of  Physical  Chemistry.  1977;66:4758. 

Welty  JR,  Wicks  CE,  Wilson  RE.  Fundamentals  of  momentum,  heat,  and  mass  transfer.  New  York 
(NY):  John  Wiley  and  Sons;  1976. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


1-D 

AP-HTPB 

ARL 

CFD 

CONP 

DFT 

HR 

MR 

TMM 

TST 

VIST 

m 

n 

TsoxJ° 

interface 


1 -dimensional 

ammonium  perchlorate-hydroxyl-terminated  polybutadiene 

US  Army  Research  Laboratory 

computational  fluid  dynamics 

constant  pressure 

density  functional  theory 

homogeneous  reactor 

microscopic  reversibility 

trial  mechanism  method 

transition  state  theory 

variational  transition  state  theory 

mass  flux  rate 

burning  rate 

temperature  at  burner  surface/condensed-phase — gas-phase 
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